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ABSTRACT 


This  report  is  the  concluding  scientific  report  to  record  the  status  and 
progress  of  Scientific  Analytical  Investigations,  the  preparation  of  Computer 
Programs,  Data  Reduction,  and  the  development  of  mathematical  and  computer 
techniques  in  support  of  Environmental  Research  and  other  various  aspects  of 
the  physical  sciences  concerning  the  upper  atmosphere. 

During  the  year  covered  by  this  report,  the  50  completed  programs  ranged 
in  complexity  and  size  from  conversion  of  programs  from  one  language  or 
computer  system  to  another,  to  analysis  and  development  of  a  large  system 
of  analytical  programs. 
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FOREWORD 


These  computer  programs  are  the  results  o£  analytical  research  performed 

for: 


The  Analysis  and  Simulation  Branch  (SUYA) 
Computation  Center  AFCRL 
Air  Force  Cambridge  Research  Laboratories 
Bedford,  Massachusetts  01730 


The  computer  programs  contained  in  this  report  may  be  obtained  from  the 
above  organization,  upon  request  by  referencing  the  appropriate  project 
number  and  problem  number  listed  at  the  end  of  each  program  description. 
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SECTION  1 


INTRODUCTION 

This  report  documents  the  efforts  under  subject  contract. 

Each  program  documented  herein  describes  what  was  done  rather  than  a  de¬ 
tailed  program  writeup. 

Upon  request,  a  complete  program  description  may  be  made  available  to 
interested  parties  (see  FOREWORD).  The  program  descriptions  are  lengthy  and 
very  detailed,  and  may  be  used  to  follow  the  complete  flow  and  operation  of 
each  program. 
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SECTION  2 


RAINFALL,  SOIL  MOISTURE  AND  TRAFFICABILITY 
IN  THE  VICINITY  OF  SAIGON 


1.  INTRODUCTION 

Trafficability,  the  ability  of  a  soil  to  permit  the  movement  of  vehicles, 
■is  a  function  of  soil  strength.  Soil  strength  is  related  to  soil  moisture,  which 
in  turn  is  related  to  many  other  variables  such  as  soil  type,  drainage  and  vege¬ 
tation.  For  the  preparat  on  of  an  accurate  estimate  of  the  effect  of  climatology 
on  soil  strength  conditions,  a  long-time  record  of  soil  moisture  measurements  is 
required.  For  most  places  in  the  world  no  such  record  exists.  In  this  study, 
equations  developed  for  estimating  soil  moisture  in  sandy  silt  of  the  type  found 
in  low  terraces  in  the  Saigon  area  are  used,  in  conjunction  with  a  22-year  record 
of  daily  rainfall  observations,  to  generate  a  soil  moisture  record.  From  knowledge 
of  soil  moisture- strength  relations  a  soil  strength  record  is  generated.  This  re¬ 
cord  is  used  to  estimate  the  probability  of  vehicle  "go"  (the  ability  to  execute 
severe  maneuvers  over  the  soil  without  becoming  immobilized)  for  given  types  of 
vehicles  on  any  given  day  of  the  year.  The  data  generated  in  this  study  are  first 
estimates  intended  for  guidance  in  planning  military  operations. 

2.  RAINFALL 

Daily  rainfall  observations  taken  at  Saigon  over  the  22-year  period  from 
1947  through  1968  were  sorted,  by  7-day  periods,  to  obtain  the  relative  frequency 
of  daily  rainfall  for  several  categories  of  rainfall  amounts.  The  frequencies  are 
plotted  In  Figure  1.  The  frequencies  for  each  7-day  period  are  based  on  154  uoser- 
vations  (7  daily  observations  multiplied  by  22  years  of  record).  The  pronounced 
seasonal  character  of  rainfall  in  Saigon  is  clearly  evident.  From  6  Jan  to  1  Apr. 
more  than  9  out  of  avery  10  days  ara  without  rain,  while  from  2S  May  through 
20  Occ.  lass  than  3  out  of  every  10  days  are  without  rain.  From  about  10  June  to 
8  Oct.  more  than  1  in.  of  rain  can  be  expected  on  l  out  of  10  days.  During  the 
22-year  period  there  were  9  days  whan  more  than  4  in.  of  rain  was  observed. 
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Figure  1,  Relative  Frequency  of  Daily  Rainfall;  *Trace,  *0.  25,  *0.  50, 
<0.  75,  <1.  00,  and  *1.  25  in.  at  Saigon 


Since  trafflcability  is  often  affected  more  by  prolonged  rainfall  than 
by  single  showers,  Figure  2  was  prepared  to  depict  the  total  rainfall  that  might  be 
expected  over  a  5-day  period.  For  example,  the  probability  of  observing  5  con¬ 
secutive  days  without  rain  is  about  85X  in  February  and  IX  from  23  May  through 
24  Oct.  During  the  22-year  period  there  were  6  five-day  periods  when  more  than 
8  in.  of  rainfall  was  observed. 

Figure  3  depicts  average  dally  rainfall  and  average  rainfall  on  a  day 
with  more  than  a  trace  of  rain.  The  curves  were  drawn  freehand .  For  example, 


Figure  3  shows  that  almost  three  times  as  much  rain  falls  on  a  rainy  day  in  June 
or  September  as  on  a  rainy  day  in  January,  and  that  lass  rainfall  occurs  near  1 


Figure  n.  Relative  Frequency  of  Totel  Rainfall  During  s  Five-Day  Perkxfc 
•Trace.  *  1. 00,  *2. 00,  eg.  00,  <4. 00,  and  *5, 00  in.  at  Saigon 
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Figure  3.  Average  Daily  Rainfall  and  Average  Rainfall  on  a  Day  With 
More  Than  a  Trace  at  Saigon 

3.  SOIL  MOISTURE  COMPUTATIONS 

Since  there  was  no  available  long-time  record  of  soil  moisture  measure¬ 
ments  taker  in  the  vicinity  of  Saigon,  daily  soil  moisture  was  computed  with  the 
use  of  equations  developed  for  estimating  moisture  from  rainfall  amounts  (Carlson 
and  Horton,  1959),  and  from  relations  derived  from  an  analogous  soil  at  site 
PD241  in  Thailand  (Kennedy  et  f)..  1967).  The  soil,  a  sandy  silt,  was  assumed 
to  be  similar  to  that  found  in  low  terraces  near  Saigon  Measurements  have  shown 
that  even  in  the  absence  of  rainfall  over  long  periods  this  soil  retains  at  least 
0.62  in.  of  moisture  in  the  0-  to  6-in.  layer  and  0.61  in  the  6-  to  12*in  layer. 
This  is  termed  the  Field-Minimum  moisture  content.  The  Field-Maximum  moisture 
content  of  the  soil,  the  highest  recurring  moisture  the  soil  can  attain,  was 
2.38  In.  in  the  0-  to  6-in.  layer  and  2. OS  in.  in  the  6-  to  12-in.  Layer.  If 
rainfall  equalled  or  exceeded  0.10  in.  in  one  day,  moisture  was  added  to  the  soil 
(accretion)  until  the  Field-Maximum  was  reached.  If  rainfall  was  not  observed,  or 
if  less  than  0.10  In  occurred,  moisture  was  lost  from  the  soil  (depletion)  until 
the  Field-Minimum  was  reeched. 
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Observations  have  shown  that  depletion  is  a  function  of  moisture  content 
in  excess  of  the  Field-Minimum.  Figure  4  depicts  sixty-order  least-squares  fits  to 


Figure  4.  Daily  Moisture  Loss  as  a  Function  of  Moisture 
Content  in  Excess  of  Field  Minimum 


measured  depletion  from  sandy  silts  for  the  0-  to  6-in.  and  6-  to  12-in  layers. 
Equations  are  provided  in  Sections  3.2  and  3.3.  Depletion  from  the  upper  layer  is 
more  rapid  than  depletion  from  the  lower  layer.  Accretion  was  assumed  to  be  a 
simple  linear  function  of  the  amount  of  daily  rainfall. 

The  method  of  generating  daily  estimates  of  soil  moisture  from  daily 
rainfall  observations  is  as  follows: 
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3.1 


Notation 


C  =  Moisture  Content  (inches),  0-  to  6-in.  Layer 
d  ■  Day 

R  =  Rainfall  (inches) 

Z  *  Available  Storage  (inches),  0-  to  6- in.  Layer 
A  *  Constant  Coefficient 

X  *  Moisture  Content  (inches)  above  Field  Minimum,  0-  to  6-in.  Layer 
M  *  Moisture  Content  (inches),  6-  to  12-in.  Layer 

Y  =  Moisture  Content  (inches)  above  Field  Minimum,  6-  to  12- in.  Layer 
S  »  Available  Storage  (inches),  6-  to  12- in.  Layer 
3.2  Computation  of  Soil  Moisture  in  the  0-  to  6- in.  Layer 

Step  1.  Assume  that  the  soil  moisture  on  31  Dec.  1946  was 
0.62  in.  in  the  0-  to  6- in.  layer,  and  0.61  in.  in  the  6-  to  12- in.  layer 
Step  2.  Compute  the  soil  moisture  on  1  Jan.  1947  using  Eq.  (1), 

Eq.  (2),  or  Eq.  (3),  depending  upon  how  much  rain  (Rd)  was  observed 
on  1  Jan.  1947  (rain  during  24-hr.  period  previous  to  observation). 


Cd+1 

-  <=d  (Z  Aixij 

if  Rd<0.10, 

(1) 

'i=l  / 

cd+l 

*  Cd  +  0.47Rd  -  0.01 

if  (Sd  +  Zd»Rd*0.l0, 

(2) 

Cd+1 

=  Cd  +  Q.75Zd  -  0.05 

if  Rd*(Sd  +  Zd), 

(3) 

where: 

Sd  -  2.05  -  Md 
Zd  *  2.38  -  Cd 
Xd  =  Cd  -  0.62 
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6 

2  A,X  -  0.030X  +  1.894X2  -  4.103X3  +  3.257X4  -  1.107X5  +  0.133X6. 
i-1 

3.3  Computation  of  Soil  Moisture  in  the  6-  to  12- in.  Layer 

Step  3.  Compute  the  soil  moisture  on  1  Jan.  1947,  using  Eq.  (4), 
Eq.  (5),  or  Eq.  (6),  depending  upon  how  much  rain  (Rj)  was  observed 
on  1  Jan.  1947  (rain  during  24-hr.  period  previous  to  observation). 


Md+1 

v1) 

if  Rd<  0.10, 

(4) 

%fl 

-  Md  +  0.22  Rd  -  0.01 

if  (Sd  +  Zd)>  Rd»0.10, 

(5) 

Md+1 

«  Md  +  0.60Sd  -  0.02 

if  Rd*(Sd  +  Zd), 

(6) 

where: 

Yd  “  1%  -  0.61 

2  Aj  Y1  =  -0.054Y  -  0.083Y2  +  3.057Y3  -  6.625Y4  +  5.069Y5  -  1.316Y6. 

3.4  Generation  of  Soil  Moisture  Record  and  Frequencies 

Step  4.  Compute  soil  moisture  for  0-  to  6-in.  and  6-  to  12-in. 
layers  on  2  Jan.  1947,  using  the  rainfall  observed  on  2  Jan.  1947 
as  Rd  and  the  Cd+i  value  computed  in  Step  2  as  Cd,  and  the  Md+^ 
value  computed  in  Step  3  as  Md. 

Step  5.  Repeat  Steps  2  and  3  until  moisture  content  values  are 
generated  for  every  day  of  the  22-year  period. 


Step  6.  For  each  7-day  period  (1  to  7  Jan.,  8  to  14  Jan.,  and 
so  forth),  find  the  relative  frequency  of  days  when  the  moisture 
content  is  in  each  of  the  following  categories  for  each  of  the  two 
layers . 


(a) 

0.70 

(j) 

(b) 

0.80 

(k) 

(c) 

0.90 

(1) 

(d) 

1.00 

(m) 

(e) 

1.10 

(n) 

(f) 

1.20 

(o) 

(8) 

1.30 

(P) 

(h) 

1.40 

(q) 

(i) 

1.50 

(r) 

1.60 

1.70 

1.80 

1.90 

2.00 

2.10 

2.20 

2.30  Omit  for  6-  to  12-in.  layer 
2.40 


4.  SOIL  MOISTURE  PROBABILITIES 

Daily  soil  moisture  estimates  computed  for  the  22-year  period  from  1947 
through  1968  were  sorted,  by  7-day  periods,  to  obtain  the  relative  frequency  of 
several  categories  of  soil  moisture  amounts.  Figure  5  shows  the  frequencies  for 
the  0-  to  6-in.  layer  and  Figure  6  shows  the  frequencies  for  the  6-  to  12-in  layer. 
Frequencies  for  each  7-day  period  are  based  on  154  estimates  of  soil  moisture. 
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Figure  5.  Relative  Frequency  of  Moisture  in  the  0-  to  6-ln.  Layer  of 
Soil;  SO.  8,  *1.  2,  *1.6,  sa.0,  «2.  2  and  *2,3in. 


2-7 


IS  3)  IS  »  IS  31  is  30  is  31  l9  30  ,9  31  |S  31  tS  30  IS  31  l8  30  |S  31 
JAN  FEB  MAR  APR  MAY  JUN  JUL  AUG  SEP  OCT  NOV  DEC 

figure  8.  Relative  Frequencies  of  Moisture  and  Rating  Cone  Index  (RCI) 
in  the  6- to  12-in.  Layer  of  Soil.  Moisture:  so.  8,  *1.2,  *1.6,  *1.8, 

*1. !),  and  *2.  0  in.  RCI:  a 50,  *77,  *120,  and  *335 

5.  SOIL  STRENGTH- VEHICLE  MOBILITY  REIATI0NS 

Traff icability  is  a  function  of  the  mass  shear  strength  of  a  critical  soil 
layer.  The  mass  shear  strength  is  expressed  in  terms  of  rating  cone  index  (RCI), 
which  is  a  measure  of  the  probable  strength  of  the  soil  under  a  moving  vehicle. 

The  depth  of  the  critical  layer  varies  with  the  weight  and  type  of  vehicle  and  soil 
profile,  but  it  is  generally  the  layer  located  6  to  12  in.  below  the  surface.  The 
RCI  may  be  compared  to  the  minimum  cone  index  required  for  50  passes  of  a  specific 
vehicle,  called  its  vehicle  cone  index  (VCI).  If  the  RCI  of  a  soil  is  higher  than 
the  VCI  of  a  particular  vehicle,  50  such  vehicles  can  be  expected  to  travel  success¬ 
fully  in  the  same  straight-line  path,  or  one  vehicle  can  be  expected  to  execute 
severe  maneuvers  without  becoming  immobilized.  A  system  has  been  developed  for 
classifying  vehicles  according  to  their  VCI  (Departments  of  the  Army  and  Air  Force, 
1968) .  A  condensation  of  the  system  for  standard  military  vehicles  is  shown  in 
Table  1. 
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Table  1.  Condensed  Classification  of  Vehicles  According  to  Vehicle  Cone  Index 
(VCI)  (contd) 


Category 

VCI  Range 

Vehicle  and  Vehicle  Typea 

2 

30-49 

Engineer  and  high-speed  tractors  with  comparatively 
wide  tracks  and  low  contact  pressures.  Examples: 

VCI  of  D7  engineer  tractor  =  40;  VCI  of  MI  14  armored 
personnel  carrier  =  37. 

3 

50-59 

Tractors  with  average  contact  pressures,  tanks  with 
comparatively  low  contact  pressures,  and  some  trailed 
vehicles  with  very  low  contact  pressures.  Example: 

VCI  of  M4B  medium  tank  -  52. 

4 

60-69 

Most  medium  tanks,  tractors  with  high  contact  pres-  1 

sures,  and  all- wheel -drive  trucks  and  trailed  vehicles 
with  low  contact  pressures.  Example:  VCI  of  M135, 
2-l/2-ton  truck  =  62. 

5 

70-79 

Most  all-wheel-drive  trucks,  a  great  number  of  trailed  . 
vehicles,  and  heavy  tanks.  Example:  VCI  of  l-l/2-ton, ; 
4x4  dump  truck  =  73. 

6 

80-99 

A  great  number  of  all-wheel-drive  and  rear-wheel- 
drive  trucks,  and  trailed  vehicles  intended  primarily  for 
highway  use.  Example:  VCI  of  1/2-ton,  4x2  pickup 
truck  =  88. 

7 

100  or 
greater 

Rear-wheel-drive  vehicles  and  others  that  generally  are 
not  expected  to  operate  off  roads,  especially  in  wet 
soils.  Example:  VCI  of  5-ton,  4x2  dump  truck  =  119. 

A  listing  of  vehicle  cone  indexes  for  military  vehicles  and  a  discussion 
of  the  procedures  for  computing  vehicle  cone  indexes,  and  a  more  detailed  discussion 
of  traff icability  terms  are  included  in  Departments  of  the  Army  and  Air  Force  <'1968). 

6.  SOIL  MOISTURE  VERSUS  SOIL  STRENGTH 

The  relation  between  soil  moisture  and  soil  strength,  in  terms  of  rating 
cone  index,  for  the  6-  to  12-in.  layer  of  soil  is  shown  in  Table  2. 

The  data  in  Table  2  were  derived  from  measurements  at  the  site  in  Thailand 
that  was  used  in  the  derivation  of  the  moisture  relations  discussed  in  Section  3, 

Soil  Moisture  Computations.  If  RCI  from  Table  2  is  substituted  for  moisture  values 
in  Figure  6  the  graph  can  be  used  to  show  the  relative  frequency  of  RCI  equal  to  or 
greater  than  a  given  value.  Thus,  during  the  rainy  period  from  the  eud  of  May  through 
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Table  2.  Relation  Between  Soil  Moisture  and  Soil  Strength  in  Terms 
of  Rating  Cone  Index  (RCI)  for  the  6-  to  12-lm  Layer  of  Soil 


Water  Content  (in.) 

6-  to  12-in.  Layer 

- 1 

RCI 

<  1.  50 

>  590 

1.50 

590 

1.50 

440 

1.60 

335 

1.65 

260 

1.70 

198 

1.75 

155 

1.  80 

120 

1.05 

95 

1.00 

77 

1.95 

62 

2.  00 

50 

2.  05 

i  40 

the  end  of  September,  moistures  of  2.0  in.  or  less  or  RCl's  of  50  or  more  occur 
80  to  90  percent  of  the  time,  and  moistures  of  1.8  in  or  less  or  RCl's  of  120 
or  more  occur  10  to  30  per  cent  of  the  time. 

7.  PROBABILITY  OF  VEHICLE  "GO" 

The  RCI  may  be  related  to  the  minimum  cone  index  requirements  of  "go" 
for  a  vehicle,  and  so  in  Figure  7  graphs  of  VCI  were  substituted  for  the  graphs 
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of  RCI  shown  in  Figure  6.  Figure  7  shows  that  during  the  rainy  season  vehicles 
with  a  VCI  of  50  will  have  an  80-  to  90-  percent  probability  of  "go"  (or  conversely, 
10-  to  20-  percent  probability  of  "no  go")  and  that  vehicles  with  a  VCI  of  120  will 
have  a  10-  to  30-  percent  probability  of  "go"  (or  conversely,  a  70-  to  90-  percent 
probability  of  "no  go"). 

8.  LIMITATIONS  OF  STUDY 

Caution  should  be  exercised  in  drawing  conclusions  on  traf ficability  in  the 
Saigon  area  from  the  graphs  presented  herein.  The  graphs  were  developed  for  only 
one  of  several  soils  in  the  area,  specifically  a  sandy  silt  found  in  low  terraces. 
Other  soils  within  the  area  include:  (1)  alluvial  and  marine  clays,  similar  to 
the  soils  in  the  lowlands  of  the  Mekong  delta;  these  clays  remain  wet  and  have  poor 

trafficability  throughout  the  rainy  season;  (2)  clays  and  sandy  clays  on  hills; 

* 

these  clays  generally  have  good  trafficability  even  during  the  rainy  season;  and 
(3)  silty  clays  in  narrow  valleys  on  terraces  and  hills;  these  clays  are  traffi¬ 
cability  problems  when  wet.  The  area  also  includes  clays  and  organic  soils  of 
mangrove  swamps  and  marshes,  which  remain  soft,  and  sands  of  river  terraces  and 
ridges,  which  remain  firm  throughout  the  year.  These  other  soils  were  not  analyzed 
because  of  limitations  of  time  and  money.  The  sandy  silt  terrace  soil  was  chosen 
for  analysis  because  it  is  one  of  the  predominant  soils  in  the  area,  and  because, 
of  all  the  predominant  soils,  it  provides  one  of  the  lowest  bearing  strengths 
(RCI)  under  wettest  conditions. 

It  should  be  recognized  that  the  charts  in  this  paper  are  based  on  rain¬ 
fall  and  therefore  are  not  applicable  to  soils  subjected  to  Irrigation.  Areas  of 
Irrigated  rice  paddles,  commonly  found  on  low  terraces,  would  normally  have  higher 
percentages  of  high  moistures  and  low  strengths  and  a  lower  probability  of  "go" 
for  a  given  vehicle  than  those  shown  on  the  charts. 


2-11 


The  probabilities  of  vehicle  "go  -  no-go"  are  applicable  only  to  soil 
conditions  on  flat  terrain.  The  slope  and  such  natural  obstacles  as  drainageways , 
scarps,  dense  vegetation,  and  microrelief  features  as  well  as  such  man-made  ob¬ 
stacles  as  railroad  embankments,  canals,  and  so  forth,  are  not  considered  in  the 
mobility  evaluations. 
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SECTION  3 


CORREIA TIONS  BETWEEN  AREAL  PRECIPITATION 
AND  GEOPOTENTLAL  HEIGHT 

This  study  was  performed  to  identify  predictors  which  might  be  used  to 
replace,  or  supplement,  persistence  in  conditional  climatologies  of  precipi¬ 
tation,  and  to  provide  at  least  partial  answers  to  the  following  questions: 

a.  Are  geopotential  height  observations  better  predictors  of  daily 
areal  precipitation  than  persistence? 

b.  Which  geopotential  height  observations  are  most  highly  correlated 

with  areal  precipitation? 

c.  How  do  predictors  of  areal  precipitation  vary  with  respect  to  the 
geographical  location  of  the  predictant  stations? 

d.  Over  how  long  a  period  of  time  are  height  observations  significantly 
correlated  with  areal  precipitation? 

To  accomplish  this  end,  computer  programs  were  written  to  correlate  geo- 
potential  heights,  geostrophic  vorticities,  and  thicknesses  between  surfaces, 
at  499  grid-points  over  one-half  of  the  Northern  Hemisphere,  extending  from 
10  degrees  E  westward  to  170  degrees  W,  at  200,  500,  and  850-mb  with  precipi¬ 
tation  observed  over  areas  in  California,  Central  United  States,  and  Eastern 
United  States  on  401  January  and  February  days. 

The  dally  total  precipitations  observed  between  0000  1ST  and  2400  1ST, 
at  10  stations  in  each  of  three  areas,  was  used  as  the  predictand.  The 
stations  are  listed  on  the  following  page  and  the  areas  are  shown  in 
Figure  1. 
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EAST 

CENTRAL 

CALIFORNIA 

1. 

Concord,  N.  H. 

Burlington,  la. 

Rfld  Bluff 

2. 

Boston,  Mass. 

DeaMoinas,  la. 

Chico 

3. 

Hartford,  Conn. 

Dubuque,  la. 

Colusa 

4. 

Albany,  N.  Y. 

Coltabla,  Mo. 

Brooks 

5. 

New  York,  N.  Y. 

Kansas  City,  Mo. 

Sacraeento 

6. 

Harrisburg,  Pa. 

8t.  Louis,  Mo. 

Stockton 

7. 

Philadelphia,  Pa. 

Springflald,  Mo. 

Modesto 

8. 

Atlantic. City,  N.  J. 

Peoria,  111. 

Merced 

9. 

Baltiaora,  Md. 

Wichita,  Kan. 

Fresno 

10. 

Washington,  D.  C. 

Tulsa,  Okla. 

Angioia  (Corcoran) 

The  data  include  all  January  and  February  days  froa  1962  through  1969. 
The  total  pradpltatlon  obaarvad  at  tha  10  stations,  la  each  area,  dur¬ 
ing  tha  period  froa  1  January  through  22  February  tor  the  yearn  1962  through 
1969 (  a  eanple  of  426  days,  was  correlated  with  the  precipitation  observed  froa 
one  to  a  lx  days  later.  The  results  vent 

OATS  1AST  CBWL  CALIFQMtU 

1  0.01  0.20  0.53 

2  -0*09  -0.00  0.20 

3  -0.03  -0.01  0.16 

4  -0.04  -0.08  0.12 

5-  -0.02  -0.03  0.11 

-0.06  0.04  0.14 


Heighcs  of  the  200,  500  and  850-tnb  surfaces  at  the  499  grid  points  shown 
in  Figure  1  were  studies.  Data  was  available  for  401  of  the  possible  424  days 
between  1  January  and  22  February  for  the  years  1962  through  1969.  Geopotential 
heights  at  each  surface,  height  difference  (thickness)  between  surfaces,  and 
geostrophlc  vorticities  were  correlated  with  precipitation.  The  vorticities  were 
simply  the  sum  of  the  heights  at  the  four  surrounding  grid  points  minus  four  times 
the  height  at  the  center  point. 

More  than  90,000  correlation  coefficients  were  computed  to  prepare  189 
maps  depicting  relationships  between  height  data  and  subsequent  precipitation. 
Examples  of  these  maps  are  shown  in  Figures  2,  3,  4,  and  5.  The  highest  corre- 
lation  coefficient  found  on  each  map  is  given  in  Table  1. 

Since  precipitation  values  used  in  this  study  are  not  normally  distributed 
in  a  statistical  sense,  the  method  proposed  by  Lund  (1970)  was  used  to  test  the 
significance  of  the  coefficients. 

The  frequency  distribution  of  the  424  dally  California  precipitation 
values  is  shown  in  Figure  6.  Thi  values  greater  than  sero  were  fitted  to  a 
normal  distribution.  The  resulting  curve, 

Y  •  0.301X-4. 378x2+12. 817X3-9.065X4+2.656X5-0.279X6  when  X  >  0.4  and  If  ■  0, 

when  X<  0.4 

where ^  is  an  estimate  of  California  precipitation  and  X  is  a  normal  number,  is 
also  shown  in  Figure  6.  Equation  above  vac  used  to  generate  8020  bogus  precipi¬ 
tation  values  (see  Figure  7).  Also,  see  Figure  8  for  the  plot  of  the  largest  of 
the  499  coefficients  obtained  with  each  of  the  20  sets  of  bogus  values. 

The  frequency  distribution  of  the  424  daily  East  precipitation  values  and 
a  least  squares  curve  fit  to  the  data, 
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Y  -  0.103+0. 173X1+1.237X13+O.171Xi4-0.251X15+O.0AlX16,  whan  2^  -0.3,  and 

Y  ■  0,  when  X  -0.3,  ^ 

a'  •'  "  ' 

where  Y  is  an  estimate  of  East  precipitation  and  is  a  normal  number  is  shown 
in  Figure  9.  See  Figure  10  for  the  highest  correlation  coefficient  obtained  from 
the  20  sets  of  499  coefficients  each  when  bogus  values  were  used*  Results  from 
Central  United  States  were  obtained  in  a  similar  manner. 

The  following  tentative  answers  were  received  from  this  study. 

1.  A  single  well-chosen  geopotential  height  observation  is  a  better  pre¬ 
dictor  of  future  daily  precipitation  than  persistence. 

2.  850-mb  heights  are  usually  better  predictors  than  observations  taken 
at  higher  altitudes. 

3.  Predictors  are  found  in  areas  expected  by  experienced  synoptic  weather 
{precasberf  rather  thim  8aM  dietsAt  un«nMpectid  idcftli^v  ^ 

4.  The  statistical  significance  of  correlations  between  height  observa¬ 
tions  and  California  precipitation  exceeds  the  5Z  level  for  all  periods 
of  less  than  five  days. 

Project  Number i  8624 

Problem  Number:  1446 
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Correlations  between  07CCL3T  geopotential  heights  (left,  .volum),  gccstrophic  vorticities  (center 
column),  height  differences  (right  column)  and  total  precipitation  observed  between  CCCOLST 
and  2UC0LST  at  10  stations  in  Eastern  United  Ltates  (darkened  qren),  on  lid  January  and 
Pctruary  days* 
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Correlations  between  C70CLST  geopotcntiel  heights  (left  column),  geos trophic  vorticitics 
(center  coluber),  height  differences  (right  column )  and  total  proeipltati.cn  observed  between 
CGCOI.OT  and  2J;GCI£T  two  days  later  (hi  to  65  hours  after  the  heights  were  observed)  at  10 
s tat inns  in  !r  sc  tern  L  ni  ted  States  (darkened  area),  on  1:01  January  and  leVnuary  days. 


Table  1.  Tr.c  hiL'h<.  ::t  correlation  coefficient  observed  on  each  oi  the  189  naps. 
The  largest  value  in  each  row  is  underlined. 


Area 
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THTCKKr-CS  ( 

T) 

200  5C0 
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2  CO 

5C0 

650 

200 

850- 

2  GO 

5G0 

“See 

650 
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-.36  .33 

—  .lih 

-.33 

.36 

.51 

.11 

.31 

.39 

Co  r» 

0 

-  .!:1  ■  .39 

.35 

-.35 

.26 

.33 
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.35 
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-.37  —61: 

-.69 
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tr 

w 
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.26 
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-.18 

.21 

Central 
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-.19  — 5U 
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-.17 
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2 
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3 
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SECTION  4 


STUDIES  OF  CLOUDS  AND  WEATHER  OVER  SEA 

I.  Introduction 

The  aim  of  these  studies  is  to  improve  the  forecasting  of  cloudiness 
and  weather  occurences  over  Southeast  Asia  (SEA)  during  the  southwest  monsoon 
enabled  by  the  supplemental  use  of  data  derived  from  satellites.  Studies 
range  from  the  development  of  simple  empirical  relationships  which  can 
currently  be  applied  in  the  preparation  of  forecasts  to  more  basic  studies 
which  are  required  for  understanding  the  weather  of  the  area. 

II.  The  Weather  Radar  Index  &  Its  Use 

A  radar  index,  RI,  was  developed  and  calculated  for  each  hour  to  provide 
a  measure  of  precipitation  activity  over  SEA.  The  index  as  calculated  in 
1969  represented  the  percental  coverage  of  radar  echoes  whose  radius  was  50 
nautical  miles  centered  at  the  four  stations  of  Tan  Son  Nhut,  Ubon,  Udorn, 
and  Plelku.  To  assess  the  accuracy  of  the  index,  a  comparison  between  computed 
areas  of  echoes  as  photographed  on  the  PPI  scope  and  computed  areas  from  radar 
echoes  (RAREPS)  was  Initiated.  This  comparison  was  first  done  using  data  from 
Miami  during  sumser  showery  conditions;  and,  later  polarcid  transparencies 
of  ground  clutter  patterns  and  the  PPI  scope  at  antenna  elevations  of  0  degrees 
rad  4  degrees  were  made  available  from  Ubon  and  Udorn  during  an  eight  day 
period  in  September  1969.  The  Miami  data  were  supplied  as  scope  photography 
tracings  which  were  ^lanimeterised  to  determine  total  echo  area  within 
a  50  nautical  mile  radius  of  the  radar.  The  SEA  data  was  projected  and  drawn 
onto  computer  generated  plots  of  the  appropriate  RAREPS.  Areas  of  the  outlined 
echoes  were  then  calculated  by  using  an  area  calculator,  which  is  an  electrical 
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device  which  counts  the  number  of  small  squares  intercepted  by  a  conducting 
pacil  as  the  operator  runs  it  over  the  areas.  From  correlation  coefficients 
between  the  measured  and  reported  indices:  0.69  for  60  cases  at  Miami, 

0.63  for  113  cases  at  Ubon,  and  0.66  for  150  c&s^s  at  Udorn,  the  regression 
lines  found  for  each  relationship  were  as  follows: 


Miami: 

RIphoto 

=  0.56 

<RIRAREP) 

+  6.96 

Ubon: 

R^photo 

*  0.46 

<RIRAREp) 

+  2.04 

Udorn: 

RIphoto 

=  0.28 

<RIRAREP^ 

+  0.64 

Thus,  on  the  average  the  RAREP  calculations  grossly  overestimate  the 
true  RI. 

In  order  to  test  for  human  variability  when  taking  RAREP  readings,  an 
analysis  of  variance  was  done  on  the  differences  between  the  RAREP  radar  index 
and  the  photography  radar  index  for  six  operators  at  Ubon.  The  results  showed 
that  one  operator  underestimated  the  radar  coverage  by  an  average  of  27.  of  the 
50  nautical  mile  circle,  while  other  operators  overestimated  by  averages  as 
great  as  87..  Using  the  F  distribution  the  averages  for  the  six  operators  are 
found  to  be  different  at  the  17.  level  of  significance.  These  differences  may 
not  be  due  solely  to  the  individual  observers  interpretations  of  observation 
but  could  arise  from  the  time  interval  between  RAREP  and  photograph  being  often 
as  much  as  15-20  minutes.  The  gain  settings  and  photographic  exposure  were 
presumed  reasonably  consistent  although  since  the  true  values  are  uncommon 
this  could  allow  for  another  degree  of  variability. 
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In  order  to  eliminate  any  significant  errors  from  Indices  computed  from 
the  RAREPS  and  since  the  true  index  varied  considerably  front  hour  to  hour,  the 
index  was  averaged  for  use  over  twelve-hourly  or  at  least  three-hourly  periods. 
Previously,  an  RI  was  emphasized  as  an  index  of  aerial  coverage  of  rain  and  not 
a  good  index  of  the  amount  of  rainfall  within  the  area.  Thus,  a  study  was  con¬ 
ducted  using  the  Tan  Son  Nhut  precipitation  record  end  Rl'e  for  that  station  over 
the  stumer  months  of  1966-68.  On  an  hourly  basis,  a  correlation  of  0*72  between 
RI  vs.  rain  occurence  was  found  using  a  sample  of  244  cases.  The  correlation  was 
0.68  for  RI  vs.  the  number  of  hours  of  rainfall  leas  than  or  equal  to  0.0  me;  and 
a  correlation  of  0«58  vs.  the  number  of  hours  of  rainfall  over  0.5  am.  When  the 
data  from  four  stations  within  50  miles  were  used,  a  correlation  of  0.82  was  found 
between  the  RI  and  hours  of  rain,  thus  substantiating  the  view  that  an  RI  rep¬ 
resents  the  frequency  or  probability  of  rain  better  than  the  amount  of  rainfall. 

From  these  findings,  the  degree  of  activity  with  regard  to  aerial  rainfall 
coverage  is  expressed  as  a  tertile  1,  2,  or  3  day.  The  part  of  the  day  consid¬ 
ered  is  from  123Q  to  2330  LT,  although  other  periods  could  have  been  used  if  de¬ 
sired.  The  average  RI's  for  this  12-hourly  period  are  determined  for  each  station, 
Tan  Son  Nhut,  Plalku,  Ubon  and  Uricra.  September's  RIs  have  been  added  for  1967  and 
1968  at  Pleiku,  Ubon,  and  Udom;  while  Ten  Son  Nhut  also  includes  June  through 
August  1966.  The  average  dally  Rl'e  fer  each  station  are  then  ranked  and  divided 
into  thirds.  The  highest  third,  TERTILE  1,  represents  active  days  while  days  in 
TERTILE  2  are  considered  average  and  days  in  TERTILE  3  are  considered  inactive. 
Since  these  classes  are  relative,  the  large  average  over-estimates  of  RI  that 
are  known  to  exist  are  unimportant. 
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Ill 


Relationship  Between  Cloudiness .  Current  and  Subsequent  Radar  Indict* 


New  correlation  coefficients  of  noon tins  cloudiness  and  coordinates  of 
the  maximum  correlation  values. 

Expansion  of  the  data  sample  to  include  Septembers  of  1967  and  1968  re* 
quircd  that  noontime  APT  pictures  for  these  months  be  prepared  and  read.  The  data 
were  then  smoothed  over  4*  X  4°  overlapping  squares  as  before,  added  to  the  accu¬ 
mulated  sample  of  1969  data  and  correlated  with  subsequent  radar  Indices  as  before. 
New  coordinates  of  the  center  of  the  4*  X  4s  squares  which  gave  the  highest  corre¬ 
lations  at  key  locations,  the  value  of  the  correlation  coefficient  and  number  of 
days  used  are  shown  in  Table  1.  This  data  is  one  of  several  variables  used  to 
yield  a  forecast  of  activity. 

IV.  Relationships  Between  Large  Scale  Pressure  Distribution  and  Subsequent 

Radar  Index 

Previous  analysis  on  this  data  was  done  using  an  average  RI  for  the  SEA 
peninsula  area.  Each  days'  1230-2330  LT  average  RI  was  an  average  of  the  combined 
RI  of  Tan  Son  Nhut,  Pleiku,  Ubon  and  Udorn.  However,  the  findings  from  day  to 
day  correlations  between  the  widely  separated  stations  were  as  low  as  +0.16,  thus 
to  possibly  secure  more  significant  results,  all  analysis  was  performed  with  re¬ 
spect  to  activity  at  a  single  station  at  a  time.  Analysis  was  done  in  detail  for 
Tan  Son  Nhut  and  to  a  lesser  extent  for  Udorn.  Also,  the  data  used  in  the  analysis 
was  expanded  to  include  September  of  1967  and  1968.  The  area  of  the  analysis 
was  also  expanded  to  cover  the  Northern  Hemisphere  in  order  to  show  planetary  wave 
patterns  if  they  truly  exist.  Other  new  date  consisted  of  QOZ  MSL  (mean  sea  level) 
grid  point  pressures  that  covered  an  area  extending  from  about  80°E  to  160SE  and 
0°  to  15*N.  These  pressures  were  read  to  0.2  millibars  from  maps  prepared  by  the 
Hong  Kong  Royal  Observatory,  and  were  then  adjusted  to  heights  at  1000  millibars 
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to  provide  an  aarial  ax tans ion  to  tha  hamispharic  1000  mb  height  analysis.  The 
MSL  data  was  evaluated  separately  to  obtain  printout  grids  which  were  then  com¬ 
bined  with  the  1000  millibar  printout  grlda. 

Generally,  1970's  results  concur  with  previous  findings.  Of  greatest  sig¬ 
nificance  was  the  feet  that  concurrently  and  for  several  days  preceding  activity 
over  the  peninsula  as  a  whole,  and  also  at  individual  stations,  heights  are  high 
over  North  Korea,  especially  at  the  200  mb  level.  Pressure  height  thicknesses 
over  chat  area  indicate  that  the  entire  air  column  is  2  1/2°  to  3s  warmer 
before  active  than  before  inactive  days.  Present  analysis,  however,  failed  to 
show  any  relationships  with  mid-latitude  troughs  at  any  level  or  abnormal  effects 
over  or  near  the  Tibetian  Peninsula.  Preceding  activity,  the  hemispheric  analysis 
yielded  an  unexplainable  picture  of  slightly  higher  heights  in  the  lower  atmos¬ 
phere  around  the  hemisphere  near  20*N  except  from  10°  to  100°E.  At  the  500  mb 
level,  heights  near  30aN  are  slightly  above  average  around  the  globe  on  active  days, 
while,  at  the  same  time  the  polar  vortex  shifts  slightly  toward  the  North  American, 
northern  Atlantic  and  northwestern  European  side.  These  departures  from  average 
are  at  most  only  2  dekameters  so  their  significance  is  uncertain  at  this  time. 

Sea  level  data,  however,  shows  a  clear  relationship  between  activity  produced  on 
the  peninsula  and  the  movement  of  areas  of  lower  pressure  westward  across  the 
Western  Pacific  Ocean  to  position  in  the  South  China  Sea  or  over  the  Philllppines. 
Figure  1  illustrates  this  westward  movement  by  the  positions  of  the  maximun  corre¬ 
lation  between  the  1000  mb.  heights  6,  5,  4,  3,  2  and  1  days  before  the  RI 
at  Pleiku.  Correlation. coefficients  at  these  points  were  about  -0.5.  For  activity 
at  Tan  Son  Nhut  the  track  of  slightly  lower  pressure  is  near  10*N,  while  for  ac¬ 
tivity  at  Udorn  it  is  at  15*N  over  the  China  Sea  and  at  20*N  at  140*E.  The  new 
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data  did  not  permit  back  tracking  of  these  minor  pressure  changes  east  of  about 
155°.  Other  researchers  have  previously  illustrated  the  westward  movement  of 
cloud  clusters  across  the  Pacific  Ocean  through  use  of  successive  days  of  photo¬ 
graphy.  This  was  further  investigated  adding  credibility  to  the  pressure  and 
cloudiness  correlations  found  at  points  east  of  the  peninsula.  Thicknesses  indi¬ 
cated  that  the  disturbed  areas,  on  the  average,  were  cool  at  low  levels  and  warm 
aloft.  As  these  surface  disturbances  move  westward,  low-level  pressures  also 
increase  over  the  South  China  coast  to  cause  easterly  winds  over  that  area.  The 
resulting  East-West  trough  is  called  the  monsoonal  trough  and  extends  from  Northern 
India  eastward  to  the  China  Sea.  Pressures  south  of  this  trough  also  fall  during 
activity,  so  the  pressure  gradients  south  of  the  trough  only  increase  slightly 
with  activity.  On  inactive  days,  the  trough  is  very  weak  or  missing  altogether. 

V.  Application  of  the  Above  Findings  to  Give  Objective  Forecasts  of  Future  Activity 

1.  RI  limits  corresponding  to  different  tertiles. 

RI  values  are  calculated  on  a  day-to-day  basis  in  order  to  verify  acti¬ 
vity.  Radar  echoes  are  plotted  on  polar  coordinate  paper  and  through  the  use 
of  an  area  calculator  (planiometer)  averages  are  computed  and  the  tertiles  are 
determined  from  Table  2  which  is  based  on  data  from  June-August  1966  plus 
June-September  1967  and  1968  for  Tan  Son  Nhut  and  June-September  1967  and  1968 
for  Pleiku,  Ubon  and  Udorn. 
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Table  1.  Correlations,  r,  between  noontime  coded  cloudiness  and  RI, 

number  of  pairs,  n,  and  coordinates  of  position  of  the  maximum 
correlation 

RI  predict  r  n  Center  4°x4° 

period  (LT)  position 


°N  °E 


Tan  Son  Nhut 

(based 

on  June  - 

Aug  66, 

June  -  Sept. 

67  and 

68  data) 

1230-1430 

same  day 

.32 

193 

11 

105 

1530-1730 

ii  it 

.36 

196 

10 

104 

1830-2030 

ii  ii 

.47 

184 

9 

113 

2130-2330 

II  II 

.48 

186 

9 

114 

1230-2330 

II  II 

.50 

191 

9 

113 

1230-2330 

Day  +  1 

.35 

182 

11 

116 

1230-2330 

Day  +  2 

.24 

188 

7 

113 

1230-2330 

Day  +  3 

.29 

179 

6 

114 

Pleiku  (Based 

on  June  -  Sept 

1967  and 

68  data) 

1230-1430  same  day 

.45 

151 

13 

105 

1530-1730  " 

.38 

152 

16 

109 

1830-2030  »  » 

.37 

167 

13 

105 

2130-2330  »  " 

.43 

163 

12 

106 

1230-2330  "  " 

.45 

173 

13 

105 

1230-2330  Day  +  1 

.42 

145 

13 

115 

1230-2330  Day  +  2 

.38 

143 

12 

115 

1230-2330  Day  +  3 

.36 

165 

15 

107 

Ubon  (based  on  June 

-  Sept 

1967  and  68) 

1230-1430  same  day 

.42 

156 

15 

108 

1530-1730  "  " 

.36 

159 

17 

108 

16 

1830-2030  "  " 

.40 

164 

15 

108 

2130-2330  "  " 

.37 

167 

14 

107 

1230-2330  "  " 

.46 

168 

15 

108 

1230-2330  Day  +  1 

.36 

168 

16 

108 

1230-2330  Day  +  2 

.44 

156 

18 

111 

1230-2330  Day  +  3 

.47 

157 

17 

112 
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Table  1.  Continued 


RI  predict 

r 

n 

Center  4°x  4° 

period  (LT) 

position 

°N  °E 


Udom  (based  on  June  -  Sept  1967  and  1968) 


1230-1430 

same  day 

.43 

154 

16 

104 

1530-1730 

It 

II 

.24 

154 

16 

106 

1830-2030 

II 

II 

.36 

163 

15 

105 

2130-2330 

II 

II 

.35 

163 

17 

104 

1230-2330 

II 

II 

.41 

163 

15 

105 

1230-2330 

Day  +  1 

.34 

141 

13 

114 

1230-2330 

Day  +  2 

.27 

134 

18 

116 

1230-2330 

Day  +  3 

.21 

124 

15 

118 

Table  2  Average  RI,  1230-2330LT  vs  tertile 


Tertile  1 

Tan  Son  Nhut 

>  17.70 

Pleiku 

>  11.75 

Ubon 

>  10.50 

Udorn 

>  V  .70 

Tertile  2 

Tertile  3 

7.01-17.69 

7.00 

5.41-11.74 

£  5.40 

4.06-10.49 

£-4.05 

5.61-14.69 

-^5.60 

Average  afternoon  tertiles  were  not  computed  unless  at  least  nine  of 
the  hourly  values  were  known. 

2.  Variables  tested 

In  order  to  derive  useful  equations  the  following  variables  were 
tried  in  a  stepwise  multiple  regression  program: 

a.  Noontime  cloudiness  averaged  over  a  4°x4°  square  at  a  key 
geographical  position. 

b.  The  MSL  pressure  gradient  at  00Z,  stations  4865Z-48914,  in  tenths 

of  a  mb. 

c.  The  OOZ  MSL  pressure  minus  1000,  at  a  key  grid  point,  or  station. 
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d.  Tertile  trend  (day-2  to  day-1) 

e.  Tertile  recurrence  (day-1)  to  forecast  day. 

f.  "Super"-trend  (day-3  to  day-2  to  day-1).  This  was  designed  to 
build  in  a  4  day  cycle. 

Tests  to  date  indicate  that  variables  a,  b  and  c  are  useful  for  calcu¬ 
lating  tertiles. 

3.  Equations  for  forecasting  tertiles 

From  the  data  period  June  -  September  1967  and  1968  the  following 
equations  appear  to  be  useful.  After  computation  of  a  tertile  number  to  the 
nearest  tenth,  a  table  is  entered  to  give  the  probabilities  of  any  tertile 
for  the  forecast  period;  these  values  can  also  be  compared  with  those  found 
if  recurrence  alone  is  used.  The  tertile  number  calculated  from  the  equa¬ 
tion  also  serves  to  give  some  indication  of  extreme  activity  or  inactivity 
if  it  is  very  low  or  very  high  respectively. 

FOR  TAN  SON  NHUT: 

Day  +  1 

Tertile  »  -0.53C  +  0.14  P  +  O.OUp  +  1.10 

where: 

C  ■  average  coded  noontime  cloudiness  1  over  4°x4°  square  centered 
at  12N,  USE 

P  ■  002  MSL  pressure  minus  1000  mb  at  8N,  124E 

4 p  ■  002  MSL  pressure  gradient  in  tenths  of  a  mb,  stations  48657  minus 
48914 

Probabilities  of  a  tertile  1,  2  or  3  day  from  the  forecast  tertile 
quantity  are  as  follows  probabilities  by  recurrence  day-1  to  day+1  are  in 
parentheses . 
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Probabilities 


Forecast 

Tertile 

Tert  3 
(inactive) 

Tert  2 
(average) 

Tert  1 
(active) 

>  2.2 

53%  (32) 

37%  (42) 

10%  (26) 

1.9  -  2.1 

18%  (29) 

46%  (42) 

36%  (29) 

4  1.8 

10%  (21) 

31%  (31) 

59%  (48) 

(3)  A  valid  average  cloudiness  requires  that  at  least  14  of  the  16  values  be 
known.  Reading  and  coding  the  clouds  is  described  on  p4,  4th  Prog.  Rept. 

Briefly:  1-0-0. 2,  2-0. 3-0. 7  and  3  *  0. 8-1.0  coverage.  Those  areas  which 
appear  to  consist  of  transparent  cirrus  are  not  counted. 

FOR  PLEIKU 

Day  +  1 
Tertile 

-  O.UP  -  0.42C  -  O.OUp  + 

1.68 

where  P 

■  OOZ  MSL  pressure  minus  1000  mb  at  11N,  126E 

C 

■  average  coded  noontime  cloudiness  over  a  4°x4° 
at  13N,  115E 

square  centered 

.\p  «  OOZ  MSL  pressure  gradient  in  tenths  of  a  mb, 
minus  48914 

stations  48657 

Probabilities  of  a  tertile  1,  2  or 

3  day  from  the  forecast  tertile  quantity 

are  as  follows;  probabilities  by  recurrence  day-l  to  day+1  are  in  parentheses. 

Probabilities 

Forecast 

Tertile 

Fertile  3 
(inactive) 

Tertile  2 
(average) 

Tertile  1 
(active) 

>  2.3 

647.  (52) 

2 71  (18) 

9%  (30) 

1.8  -  2.2 

32%  (34) 

30%  (26) 

38%  (40) 

f  1.7 

7%  (19) 

21%  (35) 

72%  (45) 

FOR  PUEIKU 
Day  +  2 

Tettile  -  -0.61C  +  0.10P  +  1.96 
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where: 


C  *  average  coded  noontime  cloudiness  over  a  4°x4°  square  centered 
at  12N,  USE 

P  *  00Z  MSL  pressure  minus  1000  mb  at  13N,  143E  (or  at  Guam) 
Probabilities  o£  a  tertile  1,  2  or  3  day  from  the  forecast  tertile 
quantity  are  as  follows;  probabilities  by  recurrence  day-1  to  day4-2  are  in 
parentheses . 

Probabilities 

Forecast  Tertile  3  Tertile  2  Tertile  1 

tertile _ (inactive) _ (average) _ (active) 

>  2.4  63%  (52)  24%  (30)  13%  f  18) 

2.0  -  2.3  27%  (32)  54%  (38)  19%  (30) 

1.9  23%  (25)  23%  (34)  54%  (41) 


FOR  P1EIKU 
Day  +  3 

Tertile  -  0.08P  -  0.41C  -  0.0l4p  +  2.04 

where  P  *  002  MSL  pressure  minus  1000  mb  at  13N,  143E  (or  at  Guam) 

C  ■  average  coded  noontime  cloudiness  over  a  4°x4°  square 
centered  at  15N  107E 

AP  »  002  MSL  pressure  gradient,  in  tenths  of  a  mb,  stations 
48657  minus  48914 

Probabilities  of  a  tertile  1,  2  or  3  day  from  the  forecast  tertile 
quantity  are  as  follows:  probabilities  by  recurrence  day-1  to  day+3  are 
in  parentheses . 

P-  >  .  Uties 

Forecast  '"'art  lie  3  Tertile  2  Tertile  1 

tertile _ ^inactive) _ 

>  2.2  59%  (48)  29%  (30)  127.  22) 

1.8  -  2.1  26%  (28)  36%  (40)  387  (32) 

<1.7  9%  (18)  41%  (35)  50%  (47) 
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4.  Evaluation  of  the  forecasting  equations. 

An  evaluation  of  the  skill  obtained  using  the  above  equations  to  fore¬ 
cast  a  tertile  1,  2  or  3  day  is  given  in  Table  3. 

Table  3.  Evaluation  of  equations  to  forecast  tertile  1,  2  or  3  days 


Station  and 
forecast  time 

Correl. 

Coef 

r 

n 

Z  correct 

Heidke  skill  score  over 
chance  over  recurrence 

Tan  Son  Nhut 

Day  +  1 

.49 

94 

52 

.28  .18 

Pleiku  Day  +  1 

.50 

102 

53 

.30  .21 

Day  +  2 

.49 

102 

57 

.37  .26 

Day  +  3 

.48 

114 

48 

.13  .08 

The  development  of  equations  for  other  stations  and  time  periods  is 
planned.  Unfortunately  these  equations  have  not  been  tested  on  the  inde¬ 
pendent  data  of  1969.  However  this  is  scheduled. 


S.  Forecast  example 

Forecast  period:  1230  •  2330  LT  next  day  at  Pleiku 

a.  From  002  surface  map  estimate  sea  level  pressure  at  UN,  126E 
to  0.2  mb.  Read  1010.0  mb,  use  10,0 

b.  From  noontime  APT  picture  read  1  square*  of  cloudiness  for 
area  11-1SH,  113-117E.  Average  to  nearest  0.1;  assume 
reading  is  1.0. 

c.  From  002  surface  map  read  sea  level  pressure  at  stations 
48657  and  48914.  Determine  difference  (48657-48914),  assume 
difference  •  0.0  mb. 

d.  insert  values  in  a,  b  and  c  in  appropriate  equation  to  read: 
Tertile  •  0.11  x  10.0  -  0.42  x  1.0  •  0.01  x  0.0  +  1.68 

e.  Compute:  Tertile  “  2.35 

f.  Read  forecast  probabilities  from  table: 

642  for  tertile  3  day 
271  for  tertile  2  day 
97,  for  tertile  1  day 
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/  HEIGHTS  VS  R I 
I30*E  SUBSEQUENT  DAYS 


Work  on  this  project  is  continuing. 
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SECTION  5 


WIND  ANALYSIS  -  UPPER  AIR 

A  large  volume  of  data  generated  by  weather  stations  throughout  the  Air 
Force  Network  is  being  gathered,  analyzed,  and  prepared  for  input.  The  data  is 
made  up  of  two  general  categories  and  falls  within  the  time  period  of  June  1967 
through  September  1968.  The  first  category  consists  of  pre-processed  cards  which 
have  been  edited  and  punched  at  Regis  College.  The  second  source  of  wind  informa¬ 
tion  is  a  collection  of  ETAC  weather  observation  tapes  which  have  been  generated 
at  a  central  facility.  Each  tape  is  a  composite  of  approximately  300  stations 
reporting  in  from  within  the  entire  network  of  weather  stations. 

The  primary  objective  of  researchers  at  AFCRL  is  to  collect  all  pertinent 
data  sources  and  establish  a  complete  and  accurate  data  set  for  each  station  on  a 
day-by-day  status.  A  complete  data  set  is  defined  as: 

1)  station  identification  block: 

a)  block,  number 

b)  altitude 

c)  latitude  &  longitude 

d)  time  in  hours 

e)  date 

2)  gradient 

3)  upper  air  data  points  (850  mb,  700  mb,  500  mb,  300  mb,  200  mb) 

A  data  point  is  composed  of  two  parameters,  wind  direction  and  velocity. 

The  final  objective  is  to  use  these  data  sets  as  input  to  an  existing 
wind  analysis  program. 

The  ETAC  tapes  will  be  used  as  a  secondary  or  back-up  source  of  data. 

If  the  pre-processed  data  set  is  incomplete,  then  ETAC  information  will  be  added 
to  the  data  set  in  order  to  establish  a  complete  record. 
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A  complete  and  exhaustive  analysis  of  the  EXAC  tapes  must  be  conducted. 

A  design  configuration  must  be  proposed  in  order  to  edit  and  select  the  variety  of 
data  samples  which  constitute  an  ETAC  tape  record.  Certain  guidelines  must  be 
established  with  respect  to  data  acceptance. 

Since  magnetic  tapes  will  be  the  primary  source  of  input,  agreement  with 
format  required  by  the  wind  analysis  program  must  be  strictly  adhered  to. 

Data  samples  from  both  sources  must  be  categorized  and  compared. 

Sorting  and  merge  techniques  must  be  utilized  in  generating  a  master  tape 
with  complete  and  accurate  data  sets.  The  work  performed  under  this  project  is 
continuing. 


Project  Number:  8624 

Problem  Number:  1694 

Researcher:  Mr.  J.  Conover 
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SECTION  6 


Windy  Acre* 

The  Boundary  Layer  Branch,  Meteorology  Laboratory  of  AFCRL,  conducted 
a  project  known  by  the  code  name  WINDY  ACRES.  This  problem  was  submitted  to  the 
Analysis  and  Simulation  Branch  (CRMCA)  under  contract  number  F19628-70-C-0029, 
problem  nmber  1764  and  project  number  8624. 

The  purpose  of  this  project  was  to  study  the  distribution  of  wind  speeds 
on  the  Windy  Acres  Tower  and  the  power  law  index. 

Project  WINDY  ACRES  resulted  in  the  creation  of  magnetic  tape  records 
during  the  summers  of  1967  and  1968  at  heights  of  0.5,  1.0,  2.0,  4.0,  8.0,  16.0, 
24.0  and  32.0  meters  above  ground  on  an  Instrumental  tower  located  about  35  miles 
northwest  of  Liberal,  Kansas.  The  data  was  recorded  for  all  eight  levels  at  1 
second  intervals  for  an  overall  total  of  40  hours  in  1967,  39  of  which  are  usable, 
and  50  hours  in  1968.  The  data  is  given  in  cm/sec  and  at  a  rate  of  3600  ,>er  hour. 

The  following  tasks  were  required  by  the  project: 

1.  To  find  the  distribution  of 

(a)  1- second  wlndspeeds  at  each  level. 

(b)  3-second  average  wlndspeeds  at  each  level. 

(c)  the  highest  3-second  average  wlndspeeds  (or  gusts)  in  1  minute 
at  each  level. 

2.  Referring  to  the  99th  percentile  of  the  1-second  wlndspeeds  at  height 
as  Vj,  find 

(a)  the  average  v  (i/j)  of  the  1-second  wlndspeeds  at  each  level 
(h^  »  0.5,  1.0,  *  •  •  32.0  meters)  for  those  seconds  when  the 

wlndspeeds  at  height  hj  is  v(j). 
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(b)  the  ratio  of  each  average  v(i/j)  to  v(j). 

(c)  to  estimate  by  method  of  least  squares ,  the  index  p  in  the 
"power  law"  for  each  key  level  (hj)  given  by 

y(t/j)  -/mV  1-1,8 

3.  Same  as  2.  except  using  the  3-second  average  windspeeds. 

4.  Referring  to  the  99th  percentile  of  the  highest  3-second  average  wind- 
speeds  in  1  minute  at  height  hj  as  g(j),  find 

(a)  the  average  v(i/J)  of  the  highest  3-second  average  windspeeds 

in  1  minute  at  each  level  (h^  -  0.5,  1.0,  *  •  •  32.0  meters)  for 
those  seconds  when  the  highest  3-second  average  windspeeds  at 
height  hj  is  g'J> 

m 

(b)  the  ratio  of  each  average  v(i/j)  to  g(j). 

(c)  to  estimate,  by  method  of  least  squares,  the  index  p  in  the 
"power  law"  for  each  key  level  (hj)  given  by 


The  derivation  of  the  index  p  by  method  of  least  squares  is  shown  below. 


GIVEN: 

LET: 


then 
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By  dropping  subscripts  for  simplicity  snd  clsrity  purposes, the  equation 


becomes 

A  -  Rp. 

taking  logs  of  both  sides, 

log  A  ■  p  log  R  or 

p  log  R-  log  A  ■  0  ; 

and  calling  E  the  error, 

E  ■  p  log  R  -  log  A. 

Then  by  squaring  E  and  expanding  the  right  hand  side, 

E^  ■  p^  (log  R)^  -  2p  log  R  log  A  +  (log  A)^ 

and  setting  d|£  ■  0  for  the  least  squares  error 
dp 

7  2 

dET_  -  2p  (log  R)  -  2  log  R  log  A  -  0 

dp 


•  • 

log  R  log  A 

u 

(log  R)‘ 


Remembering  that  the  subscripts  were  dropped,  by  adding  them  now  and 
substituting  back  the  result  is 


A  very  important  aspect  of  this  study  was  to  determine  the  value  of  p,  which 
relates  wlndspeeds  from  one  height  to  another  hslght.  It  was  particularly  desirable 
to  find  p  for  the  high  wlndspeeds. 
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The  ultimate  purpose  of  this  study  wes  to  improve  our  climatological 
knowledge  of  gusts  as  it  is  needed  to  determine  the  risks  due  to  strong  gusts  in 
structural  design  and  operation  of  equipment  and  buildings.  Studies  under  this 
project  are  continuing. 

Project  Number:  8624 

Problem  Number:  1764 

Researcher:  Mr.  Z.  Gringorten 
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SECTION  7 


MICROWAVE  ACOUSTICS 


V  “ 


Radar  and  communications  systems  require  rapid  acquisition  and  processing 
of  large  amounts  of  data.  For  a  long  while  this  was  handled  by  electromagnetic 
waves.  But  electromagnetic  devices  require  large  amounts  of  coaxial  cables.  The 
cables  are  unwieldy  as  well  as  inefficient.  What  is  needed  are  devices  which 
are  compact,  lossless,  and  easily  accessible.  The  field  of  microwave  acoustics 
is  able  to  satisfy  these  requirements,  originally  in  the  form  of  bulk  waves 
but  more  recently  by  surface  waves.  The  shift  from  bulk  devices  to  surface 
wave  devices  has  come  about  because  it  is  practically  impossible  to  vary  the 
delay  and  quite  difficult  to  manipulate  the  acoustic  energy  during  transmission 
through  a  bulk  device.  Surface  waves,  however,  have  all  the  desirable  proper¬ 
ties  of  bulk  waves,  as  well  as  being  quite  accessible  for  other  changes.  Micro- 
wave  acoustics  exhibit  the  potential  to  develop  miniature  microwave  components 
for  high-speed  signal  processing  and  avionics  systems.  This  project  is  involved 
with  the  propagation  characteristics  of  surface  waves  on  general  piezoelectric 
crystals  surfaces  in  the  presence  of  perfect  electric  and  magnetic  conductors. 

A  rectangular  co-ordinate  system  with  the  X^  axis  perpendicular  to  the 
crystal  surface  and  the  X,  axis  in  the  direction  of  propagation  is  selected. 
Various  orientations  of  the  crystal  surface  with  respect  to  the  crystal  axes 
are  to  be  considered.  In  order  to  accomplish  this,  the  crystal  axes  must  under¬ 
go  a  rotation  through  the  Euler  Angles.  The  matrix  defining  this  rotation  is 


fCc  i  Cc  '■>  Y-  J  “i  S**  f 


<  QiO  y f  (•t'oix  fa* 

-  Ct  *3  (A  /3 


where  ,  are  the  Euler  Angles.  The  differential  equations  for  mechanical 
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displacements  and  electric  potential  are  independent  of  the  surface  under  con¬ 
sideration,  only  the  values  of  the  coefficients  change  with  the  surface  orien¬ 
tation,  relative  to  the  crystal  axes.  The  tensor  quantities  of  interest,  i.e. 


the  e 


lastic  constants  (c  ,  piezoelectric  constants  («?.  . > )  and 


the  di¬ 


electric  constants  O^,  )  are  transformed  by: 

% 


n  • 

— 

n,S-Ljti- l 

\/  C  • 

j* 

''t  C  ^ 

3 

vi* 

/.  .  e 
i  c 

n,s.  t  - 1 

3 

<4  ■ 

%s  -  > 

e. 


as  t  cc 


where  the  primed  quantities  refer  to  a  rotated  co-ordinate  system  and  the  un¬ 
primed  to  when  the  crystal  surface  and  srystal  axes  systems  coincide. 

The  equations  for  the  components  U^,  U2,  U3  of  mechanical  displacement 
and  electric  ootential,CP  ,  are 


t 


>C 


[  *  <4  I’'1'’1 

e  lii  i  ‘f'jiU  -  o 

i  VV  -  o 

Considering  displacement  and  potential  independent  of  the  co-ordinate,  we 

(tm 


are  seeking  solutions  of  the  form: 

U,  --  p.  e  e‘ 

(f^ye  C 


4  -*  h 


Substituting  these  solutions  into  the  above  differential  equations  yields  a 
system  of  four  equations  in  the  unknowns/?.  ,  /y. ,  ^y.  In  order  for  a  non¬ 

trivial  solution  to  exist  the  determinant  of  the  coafficients  must  equal  zero, 
i.e,  given  a  Vs  our  problem  is  to  find  anp(  such  that 
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-  /V 


-  c 
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has  a  determinant  of  zero,  where  the  C's  are  the  elastic  constants,  the  C*  ’s 
the  piezoelectric  constants,  and  the  *s  the  dielectric  constants.  This 
involves  solving  an  eighth  order  equation  in  'A  .  For  each  value  which  has 
a  positive  real  part  we  obtain  values  of  /9  ,  »  /*>  »  A  •  Since  field  ampli¬ 

tudes  are  arbitrary  we  take/*y  ■  1  and  find  p,  ,  p*  >  pt ,  from  three  of  the 
remaining  four  homogeneous  equations. 

Assuming  the  surface  of  the  crystal  is  free  of  stress  the  mechanical 
boundary  conditions  are 


Now  the  mechanical  displacement  and  potential  may  be  expressed  as  a 
linear  combination  of  fields  associated  with  allowable  values  ofo( ;  thus  for 


**  u.'  u>  t --(t-AS/A) 

2L  A  p  C  C 

J" 
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while  in  the  interval  h^Xj^O  (using  equation  1  and  the  above  derived  boundary 
conditions)  we  get  .  »  t  +  -x  /  \ 

Cf-  A  P<  Cx'4(  *  J  ^  J  J 

J" 

Lastly  since  the  component  of  the  normal  must  be  continuous  across  »  0 
the  electric  displacement  is  given  by 

1  t  -  ■=>* 1  tPj  *.  < -•  '• *<  i 

Using  equations  4,  5,  and  6  we  get  the  following  set  of  homogeneous 

AOJ 

i’-^w  a  r><  m  r«u  Af*>pL  .  r«  pc  >'  /t  , 
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h*  r  0>>  1  ft  r  t?  4.  *->. 

/5  L‘es.  cv  J  r  &  £.'**+■  *  => 
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By  solving  Che  equation  obtained  by  setting  the  determinant  of  the  above 
system  equal  to  zero  we  obtain  the  allowable  surface  wave  velocities. 

As  a  by-product  of  this  problem  the  components  of  stress,  strain,  electric 
displacement,  electric  potential,  electric  field,  and  average  power  flow  are 
calculated  as  functions  of  wx^.  The  computer  program  has  the  option  to  print 
out  any  or  all  of  these  quantities.  The  applicable  equations  are  to  be  found 
in  Supplement  A. 


Project  Number : 

5635 

Problem  Number: 

1614 

Researcher: 

Gapt.  A.  Slobodnik  Jr..  l\  S.  A.  F 
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Electric  Displacement 
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Power  Flow 

%-U' 

The  flow  of  complex  mechanical  power  at  any  point  in  the  piezoelectric 
medium  is  given  in  component  form  as  follows: 


The  real  part  of  this  expression  represents  the  time  average  power  flow  at  a  point. 

Since  all  fields  decay  exponentially  in  the  X3  direction  there  is  no  net 
flow  of  real  power  in  this  direction.  Thus,  only  P^  and  P2  need  be  considered. 

The  components  of  the  total  time-average  power  flow  are  as  follows: 


0*3 

pj  -  1  Re[p1]dx3  *  Re|j?5] 

p5j  *  y  Rej^jdx-j  ■  Re^Pjj  j 
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where  Re  Pj  2  means  the  real  part  of  2*  T^e  final  expressions  for  the  complex 
mechanical  power  flow  P^  and  P2  are 


■p  J.  «-L  _ I 

1  _  J—  ^  ^ 


Tut  a  ^ 

rf‘::  ('« 


+  p;  l?  r  is)  J 


L'  '  T —  ■*  +  °<  -  . 

-<  -/  *-./  /j,  \1 

,  (<')  \  j.  ft  ( C'  ^  ‘'j’tJj 


t  p* 


(O  ✓  fts  *',(  *  I v 


SECTION  8 


THE  DIGITAL  SOUNDER 

A  BASIS  FOR  REAL  TIME  RECORDING  OF  IONOSPHERIC  PARAMETERS 

A  need  exists  for  the  development  of  techniques  which  can  specify,  ir.  real 
time,  the  conditions  of  the  ionospher  spanning  large  areas  of  the  globe.  There¬ 
fore,  the  techniques  of  ionosphere  specification  must  necessarily  be  able  to 
follow  rapid  changes  in  time  and  space  in  order  to  detect  the  development  of 
geophysical  disturbances  and  to  estimate  quickly  their  effects  on  various 
communication  and  surveillance  systems. 

The  systems,  such  as  radio  circuits,  radar  systems,  and  surveillance  systems, 
that  are  affected  by  the  conditions  of  the  ionosphere  operate  in  a  wide  band  of 
frequencies  ranging  from  VLF  to  VHF.  The  performance  of  the  systems  are  adversely 
affected  by  the  disturbances  causing  a  variety  of  changes  at  different  frequencies. 
A  real  time  determination  of  ionisation  distribution  would  make  it  possible  to 
compute  the  effects  of  these  disturbances  and  other  deviations  from  the  normal. 

This  knowledge  would  permit  the  restoration  of  proper  system  functioning  in 
some  cases  and  proper  interpretation  of  the  observed  phenomena  in  other  cases. 

Real  time  determination  of  ionization  does  not  necessarily  mean  real  time 
measurement;  although  it  is  an  aim  of  ionospheric  research  to  predict  the  ioni¬ 
zation  by  inference.  However,  the  direct  measurement  by  fast  methods  will  re¬ 
main  a  basic  requirement. 

1.  Limitations  of  Ionospheric  Data  Collection  Using  Conventional  Ionosondes 

a.  The  Electron  Density  Profile 

One  of  the  goals  of  ionospheric  data  collection  is  to  determine  the  dis¬ 
tribution  of  electron  density  versus  height  and  its  variation  in  time  and  space. 
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The  electron  density  &  profile  Is  basically  determined  by  a  radar  sounding  technique 
using  Ionospheric  sounders*  Narrow  RF-pulses  are  transmitted  and  propagated  up¬ 
wards,  Into  the  ionosphere.  Every  frequency  has  a  distinct  electron  density  that 
reflects  the  wave.  The  reflected  pulse  is  recorded  and  measured.  The  time  delay 
which  is  measured  is  used  to  compute  the  height  of  the  reflection  using  a  complex 
theory  and  knowledge  of  the  underlying  ionization.  Increasing  the  radio  frequen¬ 
cies  permits  penetration  deeper  into  the  ionosphere,  detecting  the  different  layers 
until  a  final  frequency  is  transmitted,  for  which  even  the  highest  occurring 
electron  density  no  longer  reflects.  The  wave  escapes  into  space.  Therefore,  this 
technique  cannot  be  used  to  gather  information  on  the  inolzation  above  the  maximum. 

b.  The  Interference  Problem  on  Ground  Sounders 

In  the  initial  years  of  sounder  operations  when  the  transmitted  pulse  power 
was  quite  lo  *,  daytime  ionograms  were  generally  poortr  than  those  taken  at  night 
due  to  absorption.  Increasing  the  transmitter  power  improved  the  situation.  How¬ 
ever,  since  the  number  of  commercial  radio  and  communications  systems  which  use 
the  same  frequency  bands  as  the  lonosondee  aa  well  ea  their  transmitted  power 
have  multiplied  over  the  past  20  years,  the  quality  of  the  ionograms  et  night  has 
deteriorated.  Due  to  the  high  level  of  interference,  whole  bands  are  useless  for 
ionospheric  sounding  at  night  and  data  hava  to  ba  interpolated  whan  possible. 

c.  The  Airborne  Sounder  end  its  Specif io  Interference  Problems 

All  ths  interfersnce  problems  that  must  bs  overcome  art  multiplied  in  an 
aircraft.  Nevertheless,  there  are  advantages  to  an  aircraft  carrying  a  sounder 
ovsr  ground  stations.  First  of  all,  sines  an  aircraft  is  mobils,  it  can  bs  used 
in  place*  where  it  le  difficult  or  Impossible  for  sounders  such  as  on  oceans,  in 
the  Arctic,  In  ths  tropics  or  in  dssarts.  Sscondly,  it  can  bs  ussd  to  msasura 


special  phenomena,  auch  as  those  occurring  during  an  adlpse,  a  nuclear  detonation, 
or  associated  with  oagnetospheric  phenomena  where  ground  sounders  can  be  installed 
only  with  great  difficulty  and  great  cost  if  at  all.  An  airborne  sounder  cannot 
replace  a  net  of  ground  stations;  however,  since  many  ionospheric  variations  are 
under  solar  control  and  move  much  faster  across  the  earth  than  the  fastest  aircraft. 

Additional  disadvantages  of  an  airborne  sounder  compared  with  a  ground  sounder 
are  displayed  schematically  in  Figure  1.  In  heavily  populated  areas,  an  aircraft 
is  more  exposed  to  the  man-made  interference  than  a  ground  sounder.  The  reason  is 
due  to  the  fact  that  the  ground  sounder  has  only  a  small  antenna  lobe,  while  the 
aircraft's  antenna  pattern  Includes  a  lobe  pointing  downward  which  is  comparable 
with  the  one  directed  upward.  The  airborne  antenna  system  captures  almost  twice 
the  energy  of  any  ground  sounder  because  of  the  ground  reflected  sky  wave, 
d.  Time  Consuming  Photo  Processing  and  Tedious  Evaluations 
The  problem  of  interference,  whether  man-made  or  natural,  is  great  in  the 
measurement  of  ionospheric  parameters.  It  is  important  to  point  out  that  technical 
limitations  in  the  data  recording  process  of  conventional  lonosondes  are  not  fav¬ 
orable  for  real  time  analysis*  The  analog  signals  are  recorded  on  film,  processed 
and  interpreted  by  specialists*  In  order  to  obtain  spatial  characteristics,  a 
net  of  stations  is  required.  The  problem  of  coordinating  various  equipments  in  the 
seme  degree  of  sensitivity  conditions  end  to  maintain  close  contact  with  data  tech¬ 
nicians  to  ensure  consistency  in  interpretation  is  quite  difficult  end  extremely 
costly  even  for  a  rather  thin  net  of  stations* 

2*  The  Digital  lonosonde  as  a  Solution  to  Existing  Problems 
a.  General 

la  order  to  overcome  some  of  the  problems  mentioned,  CRPC  (Boundary  Inter¬ 
actions  Branch)  at  Hans com  Air  Force  Base  has  begun  on  a  program  to  develop  systems 
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Chat  produce  ionoapheric  data  In  digital  form.  In  this  form,  large  amounts  of 
data  can  be  handled  on  a  computer. 

A  proposal  to  build  a  digital  ionosonde  for  the  aircraft  has  been  funded. 
This  equipment  has  been  named  Dlgisonda  128.  The  equipment  uses  phase  coding 
and  phase  coherent  integration  as  a  means  to  improve  signal  to  noise  ratio.  This 
is  combined  with  digital  data  procasalng  to  facilitate  high  spaed  operations. 

b.  Improvement  of  Signal  to  Noise  Ratio  by  Phase  Coherent  Integration 
and  Phase  Coding 

There  are  applied  detection  techniques  for  conventional  sounders  that 
use  the  energy  content  of  noise  and  signal,  but  ignoring  the  fact  that  the  fre¬ 
quency  is  known.  The  technique  of  phase  coherent  detection  utilizes  this  infor¬ 
mation.  Figure  2  shows  how  this  is  implemented  in  the  dlgieonde. 

The  upper  line  shows  a  sequence  of  transmitted  pulses,  that  appear  at  the 
receiver  output  as  the  "Direct  Signal"  and  a  sequence  of  ionospheric  returns*  Con¬ 
secutive  echoes  of  the  same  order  occur  at  the  receiver  output  with  exactly  the 
same  time  delay  and  phase  with  respect  to  the  transmitted  signal.  The  receiver 
output  is  sampled  by  gates  which  are  derived  from  a  frequency  standard  common  to 
transmitter  and  receiver.  They  also  have  a  precise  and  fixed  time  delay  with  re¬ 
spect  to  transmitted  pulses.  This  means,  that  the  samples  of  a  specific  ionos¬ 
pheric  echo  always  have  the  scsm  phase  for  one  specific  height  gets,  and  if  they 
are  summed  or  integrated,  enhance  each  other.  In  contrast,  the  phase  of  the  ran¬ 
dom  noise  will  jump  randomly  and  tend  to  cancel  out  in  the  integration,  if  enough 
samples  ere  need. 

A  range  of  128  height  Increments  have  been  choeea  to  cover  the  area  of 
lnceraat,  since  the  integration  requires  storage  and  only  a  limited  amount  of  data 
can  ba  stored.  Greater  height  ranges  than  the  "basic  Range"  can  be  inspected 
is  leas  detail  bv  increasing  ths  spacing  between  sample  gates. 
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Th«  Integration  process  is  shown  in  Figure  3.  A  separate  mesory  place  is 
provided  for  the  integration  of  the  amplitudes  of  each  of  the  128  height  rangea. 

The  atored  amplitude  information,  which  is  initially  a  "aero"  la  sequentially  fad 
into  a  temporary  auxllliary  storage  by  precisely  timed  aignals.  Next,  it  is 
converted  into  an  analog  level  and,  at  a  time  determined  by  the  proper  sample  gate, 
is  compered  with  the  receiver  output.  If  the  receiver  output  is  smaller  than  the 
storage  value,  this  value  will  be  modified  downwards, or  if  it  is  bigger,  it  will 
be  modified  upwards.  Phase  coherent  aignals  will  always  have  the  same  amplitude 
and  sign  in  their  respective  sample  gates,  changing  the  atored  information  in  the 
same  direction,  while  the  noise  with  random  positive  and  negative  amplitudes  of 
varying  magnitude  tends  to  cancel.  The  phase  coherent  signal  amplitude  will  be 
approximated  rapidly  in  its  specific  memory  place,  while  the  noise  decreasee  with 
l//aT  as  signal  statistics  can  show.  The  input  signal- to-nolse  ratio  of  S/N  will 
thus  be  Improved  to  f  S/N  after  M  Samples,  an  improvement  of  19  decibels  for 

the  case  of  80  integrations. 

Phase  coding  was  used  to  over coma  the  coherent  Interference  problem.  This 
type  of  interference  which  is  caused  by  radio  transmitters  can  sevsrly  disturb  the 
integration  process.  The  worst  case  of  Interference  which  is  one  with  exactly 
the  same  frequency  as  that  transmitted  by  the  sounder  is  Illustrated  by  Figure  4. 
Line  D  and  E  show  the  interference  pattern  between  pulses  and  the  unwanted  radio 
signal  for  two  specific  phases  of  the  transmitted  pulses. 

The  phase  switching  of  the  signal  phase  has  been  introduced  by  design  and 
this  knowledge  can  be  applied  in  the  summing  process  by  correcting  for  the  phase 
switch.  The  Idealised  example  shows  complete  cancellation  of  the  Interference  and 
enhancement  of  the  signal.  This  is  true  even  if  the  signal  is  only  a  sautll  frac¬ 
tion  of  the  interference  amplitude,  as  long  as  the  receiver  la  not  saturated  by 
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interference.  If  the  decoding  sequence  is  delayed  by  one  full  cycle,  it  is  pos¬ 
sible  to  observe  higher  order  multiples. 

c.  The  Digisonde  Format 

The  Digisonde  128  covers  a  frequency  range  fro®  25  KHZ  to  16  MHz  in 
640  steps  of  25KHz.  The  basic  frequency  increment  of  25  KHz  can  be  changed  by  pro¬ 
grams  to  50  or  100  KHz  increments.  Pulse  widths  of  50,  100,  or  150  microseconds 

are  available.  The  basic  integration  consists  of  80  samples  per  frequency  and' 

y 

height  Increment,  but  can  be  increased  in  2  multiples.  The  pulse  repetition  rate 
can  be  selected  as  50,  100  or  200  Hz.  The  fastest  possible  ionogram  between  0 
and  16  MHz  would  be  one  with  a  frequency  stepping  of  100  KHz,  an  integration  of 
80  samples  and  a  PRF  of  200  Hz.  It  would  require  64  seconds  to  be  recorded.  In 
addition  to  measuring  the  amplitude  in  each  height  gate  in  64  levels,  which  corres¬ 
ponds  to  a  1  decibel  resolution,  the  RF-phase  is  also  measured.  The  128  ampli¬ 
tude-phase  pairs  per  frequency  are  combined  with  time  and  operational  information 
and  are  recorded  on  1/2  inch  computer  tape  in  computer  compatible  format. 

d.  Online  Printout  with  Digital  Level  Presentation 

It  is  essential  that  an  immediate  picture  of  the  ionogram  is  available  at 
the  station  in  order  to  check  for  proper  gainsetting  and  frequency,  range  and  to 
allow  for  real  time  determination  of  specific  parameters.  Figure  5  shows  an  ex¬ 
ample  of  an  ionogram  produced  with  the  Digisonde  128.  Although  the  transmitted  power 
was  only  500  Watt  (it  will  eventually  be  30  KW)  and  the  antenna  extremely  poor, 
the  quality  of  the  ionogram  is  considered  good  and  the  digital  information  is  readily 
available  for  ionogram  interpretation.  A  weighting  function  is  used  in  the  con¬ 
struction  of  the  numerals  in  order  to  enhance  the  recognition  of  the  pattern  of  the 
ionospheric  trace  (the  dark  area  Increases  linearly  with  the  number).  This  basic¬ 
ally  describes  the  Digisonde  as  it  has  been  conceived  and  is  being  developed.  The 
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new  system  has  potentialities  of  routine  usage  on  both  the  ground  and  in  airplanes 
which  go  beyond  its  use  as  a  research  tool.  However,  this  would  require  a  sub¬ 
stantially  larger  effort  before  it  can  be  realized.  Only  a  few  of  the  supplements 
that  could  be  developed  are  outlined  here: 

e.  Real  Time  Availability  of  Data  for  Centralized  Computer  Processing 

An  extension  of  the  concept  of  recording  the  data  content  of  an  ionogram 
in  digital  form  is  to  connect  the  sounder  with  a  suitably  programmed  computer.  The 
computer  can  determine  such  ionospheric  parameters  as  the  lowest  observable  fre¬ 
quencies  fmin>  the  mimimum  virtual  heights  of  the  E  and  F  layers,  the  critical 
frequencies  f0E  and  foF2  and  possibly  parameters  of  the  sporadic  S-layer.  A  more 
complicated  process  is  the  separation  of  the  two  traces,  the  0-  and  x-  component 
of  the  ionogram,  by  the  computer  and  a  subsequent  true  height  analysis  which  con¬ 
verts  the  ionogram  into  an  electron  density  profile.  The  computer  may  also  be 
used  to  control  sounder  functions,  such  as  gainsetting,  number  of  necessary  inte¬ 
grations  and  the  selection  of  the  frequency  scan  limits. 

An  electron  density  profile  requires  measurements  only  on  a  few  selected 
frequencies,  the  distribution  of  which  depends  on  certain  ionospheric  conditions. 
This  selection  could  be  optimized  through  the  use  of-  a  computer  and  would  minimize 
the  operational  time  of  the  sounder.  This  fact  is  important  in  cases  where 
sounder  interference  poses  a  problem  to  other  systems. 

The  development  of  the  do. .Oder  printer  combination  has  opened  the  vay  for 
real  time  analysis  using  computers  that  are  not  necessarily  close  to  the  trans¬ 
mitter  sites.  The  total  data  content  of  an  ionogram  can  be  transformed  by  the 
decoder  through  selection  and  compression  into  an  amount  which  can  be  transmitted 
in  real  time  via  telephone  wires.  This  compressed  data  is  the  input  to  a  central 
computer  which  can  interpret  the  information  and  update  short  term  predictions. 
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A  net  of  similar  or  identical  sounders  that  are  all  connected  to  a  central  pro¬ 
cessor  would  ensure  standardization  of  the  systems.  The  continuous  nonitoring 
could  point  out  deficiencies  of  certain  stations  and  perhaps  calibration  pro¬ 
cesses  could  be  implemented  to  allow  for  proper  correlation  of  the  data. 

3*  Synchronized  Net  cf  Sounders 

a.  Application  of  Oblique  lonograms 

An  ionospheric  sounder  operating  at  one  location  can  be  received  by  a 
simultaneously  tracking  receiver  in  another  place.  The  oblique  ionogram  obtained 
permits  the  determination  of  all  propagation  parameters  on  the  specific  circuit. 
Parameters,  such  as  LUF  (lowest  usable  frequency),  MUF  (maximum  usable  frequency) 
and  the  different  modes  as  E  propagation,  1XF,  2XF  and  other  modes  can  be  lden- 

9 

tified  and  the  effects  on  transmission  of  data  can  be  derived.  The  condition  of 
the  ionosphere  in  the  mid-point  can  also  be  determined  through  careful  analysis. 

The  first  successful  attempt  to  record  an  ionogram  from  an  aircraft  was  made  at 
night  over  a  distance  of  about  2000  km.  The  systems  used  were  a  Granger  Trans¬ 
mitter  operated  at  Alta  Loma,  Texas  and  the  receiver  of  an  airborne  Granger  Sys¬ 
tem.  It  is  quite  possible  to  have  a  combination  of  ground  airborne  sounders  or 
2  airborne  sounders  monitor  the  ionosphere  over  inaccessible  territory. 

b.  The  Digisonde  and  Oblique  Propagation 

The  Digiso.ide  concept  lends  itself  immediately  for  oblique  soundings  for 
monitoring  the  ionospheric  conditions.  It  utilizes  the  exact  timing  of  the  functions 
and  the  controlled  synthesis  of  the  frequencies  using  highly  stable  standards. 

4.  Computer  Programs 

The  first  stage  of  a  major  programming  effort  to  handle  the  Digisonde  128 
data  tapes  has  been  completed.  The  existing  programs,  written  in  FORTRAN  IV 
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^•nguage  for  the  IBM  7094/DCS  at  AFCRL,  allows  tha  printout  of  tha  digital 
lonograaa  in  a  variety  of  formats.  Noise  suppressing  ■edifications  have  been 
incorporated  into  the  progran  and  the  results  show  the  feasibility  of  present¬ 
ing  ionospheric  parameters  in  this  digital  fora.  The  extraction  of  parameters 
has  been  attempted.  Minima  layer  heights  and  critical  frequenclea  versus  tine 
have  been  successfully  extracted  from  suitable  actual  data.  The  reformatting 
of  reduced  parameters  into  a  format  usable  by  the  Digicoder-Xerox  Printer  com¬ 
bination  has  been  accomplished.  This  permits  the  use  of  the  pattern  enhancing 
scheme  in  the  presentation  of  these  data.  Figures  6  through  12  are  examples 
of  the  printed  outputs  from  the  Digicoder-Xerox  Printer  combination  and  the 
IBM  7094/DCS. 
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AIRBORNE  SOUNDER 


PRINCIPLE  OF  PHASECOHERENT  GATING 
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Minimum  Heights  of  E"  and  F- Layer  versus  Time 
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SECTION  9 


ANALYSIS  OF  RADIO  PROPAGATION  CONDITIONS  OF  THE  IONOSPHERE 

Radio  propagation  conditions  of  the  ionosphere  are  being  recorded  con¬ 
tinuously.  The  primary  source  of  research  data  is  an  analog  magnetic  tape.  The 
purpose  of  this  particular  research  project  is  to  verify  the  data  samples  recorded 
on  the  analog  tape. 

In  order  to  accomplish  this  task,  the  analog  tape  must  be  converted 
to  a  digital  tape.  The  conversion  of  analog  to  digital  was  processed  at  the  Data 
Reduction  Facility  at  Northeastern  University.  The  format  of  the  digital  tape 
is  classified  as  "DCS-7094  type  1  format".  The  actual  data  tape  used  in  the  analysis 
was  labeled  "DOPLER"  and  it  contained  "Ionospheric  Doppler  Shift"  data  recorded  at 
the  Canadian  range  on  23rd  of  November  1968.  The  sampling  rate  is  5  KHz  (5  kilo¬ 
cycles/sec).  A  tape  record  consists  of  a  time  word  and  twenty  four  (24)  data 
values.  The  data  values  represent  results  from  two  channels,  5  and  4.  The  order 
in  which  these  data  points  were  written  on  the  tape  is  the  staggered  mode;  5,  4, 

5,  4,  etc...  12  values  of  each.  It  should  be  pointed  out  that  one  of  the  program 
requirements  was  that  as  many  as  four  (4)  channels  may  be  recorded  on  future  tapes. 
The  data  processing  program  (ATA1)  can  handle  1,  2,  3,  or  4  channels  as  generated 
per  tape  record. 

"DCS-7094  type  1  format"  will  now  be  described.  The  data  is  set  up  in 
logical  records,  each  consisting  of  two  control  words  and  a  given  number  of  frames 
of  data.  The  first  control  word  is  the  number  of  data  words  in  a  frame,  including 
the  time  word.  The  second  control  word  is  the  number  of  frames  in  the  logical 
record.  These  two  control  words  are  fixed  point  or  integer  numbers.  The  "end 
of  data"  consists  of  a  record  of  only  three  (3)  words;  "1,  1,  0".  It  is  impossible 
to  have  a  frame  of  data  with  only  one  data  word.  Such  a  condition  uniquely  iden¬ 
tifies  the  "end  of  data". 
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The  "time"  word  is  a  36  bit  word  and  it  is  the  first  word  of  each  frame. 


It  is  constructed  in  the  following  manner: 
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*SYN  (13  bit)  denotes  "synchronization"  (  1  =  not  sync,  0  *  in  sync) 

Another  program  requirement  was  that  unpacking  the  time  word  to  obtain 
individual  values  of  hours,  minutes,  seconds,  and  miliseconds. 

A  data  processing  program  (ATA1)  was  written  in  FORTRAN  IV  to  run  on 
the  Direct  Couple  System  at  AFCRL.  This  program  reads  the  digital  tape,  unpacks 
the  time  word,  searches  for  the  proper  time  segment,  converts  raw  values  ('voltage) 
to  amplitude,  flogs  values  out  of  range,  prints  converted  values  with  time,  writes 
output  record,  and  plots  amplitude  points  graphically  on  line. 

The  tape  is  read  by  means  of  a  READ  TAPE  statement  as  outlined  in  CEIR 
7090  FORTRAN  COMPATIBLE  FLOATING  POINT  memo  -  R.  Mclnemey.  Also,  the  "eid  of  data" 
and  EOF  are  checked. 

By  using  the  "MOD"  function  and  "powers  of  2",  time  information  is 
stripped  from  the  time  word. 

Input  to  the  program  requires  a  start  and  stop  time.  If  the  time  of  a 
particular  record  does  not  fall  within  the  time  interval,  the  program  either 
keeps  on  searching  for  a  starting  point  or  terminates. 

The  conversion  technique  is  based  on  input  parameters:  for  each  channel, 
the  following  Information  is  required;  upper  and  lower  tolerances,  high  and  low 
ranges.  The  raw  data  value  is  converted  by  the  equation: 

V  -  x  (rgk)  +  rgl 
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where: 


V  =  converted  value 

r  ■  raw  value  to  be  converted 

tl  ■  lower  tolerance 

rgl  =  low  range 

tk  ■  (upper  -  lower)  tolerance 

rgk  =  (high  -  low)  range 

The  converted  value  is  tested  to  see  i£  it  falls  within  the  range  interval 
specified  for  each  particular  channel.  If  it  does  not  satisfy  the  limits,  the 
value  is  flogged  at  printout. 

The  on  line  plot  scheme  is  called  SCETCH.  This  display  prints  a  two 
coordinate  graph.  The  y  coordinate  is  amplitude  value  and  the  x  coordinate  is 
the  sequential  order  of  signals. 

The  scope  of  this  task  is  to  show  that  the  analog  and  digital  data  tapes 
are  compatible.  Once  this  has  been  established,  the  digital  data  (from  ATA1  out- 
put  tape)  will  be  subjected  to  power  spectrum  and  cross  spectrum  analysis. 


Project  Number: 
Problem  Number; 
Researcher : 


5631 

1689 

Dr.  c.  Rush 
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SECTION  10 


INVESTIGATION  OF  VARIATIONS  OF  STRATOSPHERIC  AEROSOL 
THROUGH  PHOTOGRAPHIC  MEASUREMENTS 

In  the  course  of  the  atmospheric  optics  research  by  CRCA,  photographic 
twilight  measurements  for  investigations  of  variations  of  stratospheric  aerosol 
have  been  made  at  several  stations.  Several  films  each  with  about  60  exposures 
(simultaneously  in  two  wavelength  ranges)  had  to  be  evaluated.  Initially  these 
calculations  were  done  by  hand  but  it  was  found  that  this  process  was  much  too 
tedious  and  to  ease  the  workload  a  computer  program  was  written  to  perform  all 
the  necessary  calculations. 

The  program  was  designed  to  compute  the  solar  elevation  from  local  time 
and  to  convert  film  densities  to  intensities,  taking  care  of  corrections  for 
neutral  filter  steps,  aperture  setting  and  exposure  duration. 

Since  the  input  to  the  program  was  given  in  terms  of  the  transmission 
percent,  it  was  necessary  first  of  all  to  convert  from  this  film  transmission 
percentage  to  film  density.  This  is  done  by  means  of  the  graph  shown  in  Fig.  1, 
which  gives  the  relationship  between  transmission  percent  and  film  density  for 
this  particular  film  12C .  Figure  1  also  shows  how  arrays  are  set  up  giving  corre¬ 
sponding  values  for  transmission  percent  and  film  density  for  the  U  and  G  filter 
ranges.  Thus,  inserting  these  arrays  into  the  program  we  can  eas ity  find  the 
density  of  the  film  given  the  transmission  percent  through  a  simple  interpolation 
scheme. 

Since  the  area  of  concern  was  the  film  intensity,  this  final  conversion 
was  necessary  and  the  equation  is  as  follows: 

0o  •  i5o  +•  k  -  log  (T/5) 

where  Co  is  the  weighted  film  density  mean,  k  is  a  constant  depending  on  the 


aperture  setting,  and  T  Is  the  exposure  duration. 

More  recently,  a  plotting  routine  «as  added  to  the  program  for  a 
detailed  study  of  these  films.  The  plots  con.lsted  of  intensity  and  color 
versus  solar  elevation.  A  sample  plot  can  be  seen  in  Fig.  2. 
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SECTION  11 


WAVE  PROPOGATION  IN  HOT  MAGNETAPIASMAS 

Recent  work  at  AFCRL  has  been  concerned  with  deriving  dispersion 
relations  for  electromagnetic  waves  propogating  in  plasmas  immersed  in  a  magnetic 
field.  These  dispersion  relations  can  be  used  to  predict  how  waves  progogating 
in  a  plasma  will  attenuate  and  distort  as  a  function  of  the  parameters  charac¬ 
terizing  the  medium.  These  investigations  are  of  use  in  the  design  of  any  radio 
communications  system  in  which  waves  must  propogate  through  an  intervening  plasma 
medium,  such  as  in  reentry  communications  systems  and  in  ionospheric  propogation 
work.  Conventional  dispersion  relations,  such  as  the  appletos-Hartree  equations, 
neglect  a  number  of  important  properties  of  plasmas  that  can  have  an  effect  on 
wave  propagation.  The  work  at  AFCRL  is  concerned  with  deriving  dispersion  rela¬ 
tions  that  are  much  more  accurate  and  that  are  valid  over  a  much  wider  range  of 
plasma  parameters  than  previous  theoretical  formulas. 

The  dispersion  relations  are  obtained  by  solving  the  Bal-tzmann  kinetic 
equation  using  a  variety  of  Fourier  transform  techniques  ar.d  perturbation  methods. 
Graphs  are  available  of  the  propogation  constants  plotted  over  a  wide  range  of 
normalized  plasma  parameters.  The  cases  of  propogation  along  and  across  the  DC 
magnetic  field  lines  are  investigated  in  detail.  The  graphs  show  the  effect  of 
velocity  dependent  electron-neutral  collisions  and  Coulomb  collisions  on  Che 
propogation  constants  of  the  wave.  Four  graphs  are  given  at  the  end  as  a  sample 
of  the  output  generated  in  support  of  the  research  on  this  project. 

The  interaction  of  electromagnetic  waves  with  an  infinite,  homogeneous 
plasma  is  completely  described  by  Maxwell's  equations  and  the  constitutive  relations 
For  such  a  medium,  Maxwell's  equations  may  be  written  ir  the  form 

VX  E  «  *  B  (1) 

VX  H  •  C  (2) 

Aa  '  - 
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where  C  =  total  current  density.  The  constitutive  relations  express  C  in  terms  of 
D  and  £.  For  a.c.  fields,  it  is  possible  to  obtain  two  separate,  yet  equivalent, 
models  of  the  plasma  by  considering  the  medium  to  possess  either  dielectric  or 
conducting  properties.  For  the  dielectric  model,  the  electrons  and  positive  ions 
are  considered  to  form  an  induced  dipole  moment  as  the  result  of  interacting  with 
an  externally  imposed  electric  field.  The  motion  of  the  electrons  then  gives  rise 
to  a  polarization  vector  P  =  -Her  where  N  is  the  electron  density,  e  is  the  elec¬ 
tronic  charge,  and  where  the  electron  displacement  r  is  determined  from  an  equation 
of  motion  of  the  electron.  In  the  following,  a  subscript  will  denote  the  Cartesian 
component  of  a  vector.  Since  =  Pi  +  g0E^,  the  total  current  density  is 

equal  to  the  electric  displacement  current  density  .  This  model  yields  the 
constitutive  relations. 


o  «i 


Pi  + 


oEi  = 


For  the  electrical  conductor  model  of  a  plasma,  the  electrons  and  ions 
are  considered  to  be  free  charges  that  acquire  a  net  drift  velocity  under  the 
influence  of  externally  imposed  fields.  Here,  the  total  current  density  Ci  is 
the  sum  of  a  conduction  current  density  and  a  free  space  displacement 

current  density  0%  >  so  that  +  0E^.  The  constitutive  relations  are 

completed  by  the  relation:  B^  *  o^i >  since  the  diamagnetic  effect  of  a  plasma 
at  ordinary  densities  is  negligible.  When  the  ion  motion  is  neglected,  the 
conduction  current  density  is  given  by  *  -Nevj  where  the  electron  velocity  vj_  is 
determined  from  the  equation  of  motion  for  the  electron.  Hence,  the  conduction 
current  density  is  exactly  equivalent  to  the  polarization  current  density  P^  in 
the.  dielectric  model.  By  equating  the  total  current  density  in  the  dielectric  model 
to  the  total  current  density  in  the  conductor  model,  the  following  expression  is 


obtained  for  the  dielectric  constant,  K; 


K  -  1  + 


iuj  € 


(3) 
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where  K  ■  £/£c  <r  “  electrical  conductivity  and  the  electric  field  varies  in  time 
as  exp  +  ( jwt) . 

In  the  following  presentation,  there  will  be  a  number  of  assumptions  made 
about  the  interaction  of  the  electromagnetic  waves  with  the  plasma  medium.  First, 
it  will  be  assumed  that  the  plasma  is  electrically  neutral  in  the  absence  of  any 
external  perturbations.  Also,  the  plasma  will  be  assumed  to  be  weakly  ionized,  so 
that  electron  collisions  with  the  neutral  particles  will  dominate  over  electron- 
electron  and  electron-ion  collisions.  The  signal  strength  of  the  electromagnetic 
fields  will  be  low  in  order  that  the  kinetic  equations  may  be  linearized.  For 
low  signal  strengths,  nonlinear  effects  such  as  harmonic  generation  and  self-modu¬ 
lation  are  negligible.  Since  the  ions  are  much  heavier  than  the  electrons,  they 
respond  more  slowly  to  an  alternating  electric  field.  It  will  be  assumed  that  the 
signal  frequencies  are  sufficiently  high  so  that  the  ion  contribution  to  the  elec¬ 
trical  conductivity  is  negligible. 

A  fundamental  length  characterising  a  plasma  is  the  Debye  length,  /D, 


where  k  ■  Boltzmann's  constant 
and  Tg  ■  electron  temperature. 

In  considering  the  influence  of  many  particles  upon  a  test  particle  introduced 
into  the  plasma,  it  is  found  that  there  is  a  limit  to  the  number  of  particles  that 
have  to  be  considered  in  the  interaction.  Charged  particles  close  to  the  test 
particle  act  as  a  screen  between  the  test  particle  and  more  distant  charged  particles. 
The  Debye  length  is  the  length  parameter  that  characterizes  this  transition.  It 
will  be  assumed  that  the  electromagnetic  wavelength  is  much  larger  than  the  Debye 
length  p .  This  is  necessary  in  order  to  consider  the  perturbing  electromagnetic 
fields  as  macroscopic  fields.  Finally,  the  plasma  temperature  will  be  considered 
to  be  sufficiently  low  so  that  the  particle  motion  is  nonrelativist ic 
k  Te  «  me2 
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where 


m  =  electron  mass. 

It  will  be  shown  that  first  order  thermal  corrections  to  a  cold  plasma  are  propor- 

9 

tional  to  n“,  where 

0*  =  kl^  / mc^  and 
n  =  index  of  refraction 

Hence,  for  nonrelativistic  plasmas,  thermal  effects  are  important  only  when  the 
index  of  refraction  is  large. 


Special  Features  of  Hot  Plasma  Theory 

There  are  several  different  approaches  that  may  be  taken  in  describing 
the  dielectric  and/or  conduction  properties  of  a  plasma.  The  most  appropriate 
method  for  a  particular  problem  depends  upon  the  nature  of  the  information  desired. 
The  simplest  model  of  a  plasma  is  based  upon  orbit  theory.  According  to  this 
approach,  it  is  assumed  that  the  individual  particles  comprising  the  plasma  do  not 
interact  with  one  another  and  that  they  do  not  possess  a  random  thermal  motion. 
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'8) 


Substituting  Eq.  (7)  into  (3)  yields 


where  Wjj  =  (Ne  /m  Q)  is  the  electron  plasma  frequency,  and 

X  X 

Although  orbit  theory  is  not  capable  of  rigorously  taking  into  account  the  effect 
of  interparticle  collisions,  an  average  collision  frequency  v  is  often  introduced 
into  the  analysis  in  a  heuristic  fashion  by  adding  a  damping  term  m  v  to  the  left- 
hand  side  of  the  equation  of  motion  (5). 

The  effect  of  the  thermal  motion  of  the  particles  cannot  be  adequately 

taken  into  account  within  the  framework  of  orbit  theory.  The  most  direct  means 

of  including  the  effects  of  thermal  motion  is  to  retain  the  pressure  gradient 

terms  that  appear  in  the  momentum  equations  for  each  plasma  constituent.  If  the 

wave  motion  is  assumed  to  be  adiabatic,  then  the  continuity  equation  and  the 

momentum  equation  for  each  constituent,  together  with  the  electrodynamic  equations, 

form  a  determined  set  of  equations.  By  linearizing  these  equations  and  looking 

for  solutions  that  vary  in  space  and  time  as  expj (wt-k«r) ,  where  k  is  the  wave 

•** 

number,  a  dispersion  relation  between  w  and  k  may  be  obtained.  In  contrast  to  the 
results  of  orbit  theory,  the  linearized  moment  equations  yield  a  dispersion  rela¬ 
tion  that  includes  first  order  temperature  terms.  For  transverse  electromagnetic 
waves  the  thermal  corrections  are  negligible  except  when  a  d.c.  magnetic  field  is 
present  and  the  signal  frequency  is  close  to  the  electron  cyclotron  frequency.  On 
the  other  hand,  for  longitudinal  waves  (in  which  the  direction  of  propagation  is 
parallel  to  the  electric  vector)  the  effects  of  finite  temperature  are  important 
because  they  provide  a  mechanism  for  the  propagation  of  the  longitudinal  oscillations. 

According  to  the  orbit  theory  description  of  longitudinal  oscillations  in 
a  cold  plasma,  if  a  group  of  electrons  is  slightly  displaced  from  its  equilibrium 
position,  a  space-charge  field  will  immediately  develop  tending  to  drive  the 
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electrons  back  to  their  equilibrium  position.  The  momentum  they  gain  will  cause 
them  to  over- shoot  their  equilibrium  position  and  a  space  charge  field  will  again 
tend  to  restore  them  to  their  equilibrium  position.  In  a  cold  plasma,  this  os¬ 
cillation  about  their  equilibrium  position  will  occur  at  the  plasma  frequency, 
w., .  However,  if  the  electrons  are  displaced  from  their  equilibrium  position  and 
they  possess  thermal  energy,  a  pressure  gradient  will  develop,  tending  to  ac¬ 
celerate  the  electrons  in  the  same  direction  as  the  space  charge  electric  field. 

The  frequency  of  these  oscillations  in  a  warm  plasma  is  slightly  greater  than  the 
plasma  frequency,  and  the  oscillations  now  can  propagate  and  carry  energy.  These 
oscillations  are  maintained  primarily  by  the  velocity  modulation  of  a  small  group 
of  electrons  that  are  moving  rapidly  in  the  direction  of  propagation  of  the  wave. 

This  dispersion  relation  for  the  electron  oscillations  is  given  by 

w2  "  "N  2  +  k2vTH  2 

where 

k  *  wave  number 

VTH2  *  3kTe^m  • 

When  a  d.c.  magnetic  field  is  present,  then  it  is  found  that  this  longitudinal 
wave  can  propagate  along  the  magnetic  field  lines.  If  the  linearized  moment  equations 
arc  used  to  study  the  propagation  of  longitudinal  electron  waves  across  the  magnetic 
field  lines,  then  the  dispersion  relation  takes  the  form  exhibited  in  Figure  1. 

Here,  w^  ■  eB/m  Is  the  electron  cyclotron  frequency.  In  this  figure,  the  frequency 
is  normalized  to  the  electron  cyclotron  frequency  and  the  abscissa,  kvjg/wH,  is  the 
ratio  of  the  mean  radius  of  gyration  to  the  wavelength.  The  analysis  based  upon 
the  linearized  moment  equations  is  not  valid  if  the  radius  of  gyration  of  the 
electrons  is  sufficiently  large  so  that  in  a  single  gyration  an  electron  samples  a 
significant  portion  of  the  electric  field.  Hence,  tha  dispersion  curves  in  Figure  1 
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are  valid  only  for  long  wavelengths.  A  more  detailed  kinetic  theory  approach  is 
necessary  in  order  to  ascertain  how  these  curves  must  be  modified  for  larger  wave- 
numbers  . 

For  certain  types  of  collision- 
dominated  plasmas  in  a  magnetic  field,  a 
conducting  fluid  model  is  often  simpler 
and  more  useful  than  the  model  based 
upon  the  moment  equations  for  each  plasma 
constituent.  For  many  cosmic  and  astro- 
physical  plasmas,  the  time  and  leng.ii 
scales  are  so  large  and  the  frequencies 
are  so  low  that  the  motion  of  the  elec¬ 
trons  relative  to  the  ions  is  unimportant. 
The  important  question  concerns  the  motion 
of  an  element  of  conducting  fluid  viewed 
as  a  whole.  Thus,  there  is  less  infor¬ 
mation  contained  in  the  conducting  fluid 
model  of  a  plasma  as  opposed  to  the  single 
constituent  moment  equations.  The  equations 
for  the  conducting  fluid  model  may  be  de¬ 
rived  from  the  single  constituent  moment 
equations  by  writing  these  equations  in 
terms  of  a  mean  velocity  for  the  mixture  and  summing  over  all  the  constituents  of 
the  gas. 

The  most  general  approach  to  wave  propagation  in  a  plasma  includes  the  use 
of  kinetic  equations  for  each  plasma  constituent.  Since  only  relatively  high  fre- 


Figure  1.  Electi  on  Wave  Propagating 
Across  a  Magnetic  Field  Based  on  the 
Montent  liquations,  Higher  order  tem¬ 
perature  effects  cause  change  in  disper¬ 
sion  curves  1 1  ^  +  *' 
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quencies  are  considered  in  this  analysis,  only  the  kinetic  equation  for  electrons 
is  of  interest.  The  kinetic  equation  describes  the  behavior  of  the  velocity  dis¬ 
tribution  function,  f.  This  distribution  function  is  defined  as  the  density  of 
points  in  the  s’x  dimensional  phase  space  made  up  of  the  three  components  of  the 
particle's  position  plus  the  three  components  of  the  particle's  velocity.  The 
kinetic  equation  approach  must  be  used  to  describe  wave  propagation  in  a  plasma 
when  information  is  desired  about  phenomena  that  depend  explicitly  on  the  velocity 
spread  of  the  particles. 

There  are  several  important  effects  which  occur  in  a  hot  magnetoplasma 
that  are  absent  in  a  cold  plasma.  One  effect  is  transconductance,  in  which  the 
conductivity  at  a  given  point  in  the  plasma  is  not  entirely  due  to  the  drift 
velocity  of  the  electrons  acquired  by  the  presence  of  the  electric  field  at  that 
point,  but  may  depend  upon  the  effect  of  the  electric  field  at  neighboring  points. 
This  nonlocal  effect  is  due  to  the  fact  that  the  electrons,  because  of  their  ther¬ 
mal  motion,  may  move  during  one  period  of  oscillation  of  the  electric  field.  This 
nonlocal,  hot  magnetoplasma  effect  may  be  expressed  mathematically  as 


t) 


JJ  <f  ij (r  -r\  t-t')  Ej(r',  t  )  •  d3r'dt' 


ao> 


Another  important  effect  which  occurs  in  a  hot  plasma  involves  collision¬ 
less  wave  damping.  In  an  isotropic  plasma,  longitudinal  vavea  may  be  damped  when 
the  wavelength  becomes  of  the  order  of  a  Debye  length.  This  is  termed  Landau 
damping  and  is  the  result  of  energy  transfer  between  the  longitudinal  wave  and  those 
particles  travelling  near  the  wave  phase  velocity.  This  energy  transfer  (s  n  con¬ 
sequence  of  the  fact  that  particles  travelling  slightly  slower  chan  the  wave  phase 
velocity  will  gain  energy  from  the  wave  if  they  are  initially  randomly  phased  with 
respect  to  it.  Another  type  of  collisionless  damping  is  termed  cyclotron  damping. 
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and  occurs  when  a  right-hand,  circularly-polarized  wave  is  travelling  along  the 
magnetic  field  lines  at  a  signal  frequency  close  to  the  electron  cyclotron  fre¬ 
quency.  If  an  electron  is  gyrating  about  the  magnetic  field  lines  at  the  cyclo¬ 
tron  frequency  ,  and  is  travelling  along  the  magnetic  field  lines  with  a 
velocity  +  'v)^,  then  in  the  electron's  frame  of  reference  the  frequency  of  the 
wave  is  (w  f^llk),  where  (v>nk  represents  the  Doppler  sh<ft  due  to  the 
relative  motion  of  the  electron  and  the  wave.  At  the  resonant  electron  velocity 
(v>H  equals  +  (v  -u^)/k,  the  electrons  rapidly  gain  energy  at  the  expense  of 
the  electromagnetic  wave.  This  damping  results  in  a  finite  maximum  wave  number 
for  propagation  of  the  right-hand  wave  at  the  cycletron  frequency. 

A  third  important  effect  which  is  predicted  by  the  kinetic  theory  approach 
involves  a  modification  of  the  dispersion  curves  shown  in  Figure  1  for  propagation 
across  the  magnetic  field.  This  modification  is  due  to  resonant  effects  that  may 
occur  when  the  Larmor  radius  (mvT/wH)  is  an  integral  number  of  electromagnetic 
wavelengths.  It  will  be  seen  that  pass  bands  can  exist  for  longitudinal  waves 
propagating  across  the  magnetic  field  close  to  all  the  harmonics  of  the  cyclotron 
frequency. 

Work  under  this  project  is  continuing. 


Project  Number:  4642 

Problem  Number:  1550 

Researcher:  Mr.  R.  Papa 
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Graph  A  shows: 


NU+  = 


(1-Y) 


3=7^) 


vs  y  for  y *  .4  to  1.1 


Graph  B  shows: 


RHO- 


f  + 


o*ii  !«•*!* 


VS 


for 


*  .4  to  1.1 


Graph  C  shows: 


RE(N)  «  Re  (Nq) 

vs  X  for  X  ■  .4  to  1.0 


Graph  D  shows: 

1M(N)  -  In  (Nq) 

vs  X  for  X  »  .4  to  1.0 


The  parameters  characterising  the  nedlun  are  listed  in  the  upper  left 
corner  of  each  graph. 
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RHO  -  # 

-3*«.00  -29.00  -2H.00  *19.00  -IN. 00  -9.00'  -H.or  1.00  6.00  11.00 


X  •  I .  ‘.00003 

20  •  .*00000 

22  •  -.100000 

VEE  ■  1.000000 

VCl  *  .500000 

DELTA*  .000001 
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SECTION  12 


SOME  MATHEMATICAL  PROPERTIES  OF  AN 
INTEGRAL  EQUATION  FOR  THREE- BODY  SCATTERING 

In  the  past  ten  years  the  quantum  mechanical  three-body  problem  has  re¬ 
ceived  a  lot  of  attention  in  the  literature.  The  main  advance  in  treating 
this  problem  has  been  the  notion  that  the  wavefunction,  or  the  T-matrix,  may 
be  decomposed  into  a  sum  of  three  parts.  The  reason  for  making  this  decompo¬ 
sition  is  the  hope  of  obtaining  equations  which  are  tractable,  whereas  a 
straightforward  approach  leads  to  many  difficulties. 

The  two  recent  methods  for  the  three-body  problem  are  given  in  the  work 
of  Eyges*  and  Jasperse^*^,  and  the  work  of  Faddeev^  and  others-*’**.  Eyges  and 
Jasperse  decompose  the  wavefunction  directly,  and  Faddeev  and  others  decompose 
the  T-matrix  directly. 

Once  the  decomposition  is  made,  two  questions  arise  immediately.  They 
are,  (1)  what  are  the  mathematical  properties  of  the  resulting  equation?, 
and  (2)  is  the  equation  tractable  using  modern  computational  methods?  In  this 
report,  a  study  of  some  mathematical  properties  of  Eqn.  (16)  of  reference  3 
is  presented.  This  equation  describes  the  simplest  model  of  neutron-deuteron 
scattering,  which  still  retains  the  basic  features  of  a  more  complicated 
nuclear  fission  process. 

In  this  first  phase  of  the  investigation  of  Eqn.  16,  we  have  studies  the 
resulting  integral  operator  to  determine  whether  or  not  it  is  compact  as  an 
operator  on  a  certain  Banach  space.  Knowing  that  the  operator  is  compact  would 
allow  us  to  apply  the  classical  Fredholm  Theory.  This  theory  enables  us  to 
say  something  about  the  existence  and  uniqueness  of  solutions  to  functional 
equations.  Without  an  answer  to  this  existence  and  uniqueness  problem,  the 
results  of  any  numerical  attempt  to  solve  eqn.  16  would  be  at  best  difficult 
to  interpret. 


12-1 


The  results  to  be  proved  in  this  report  are  summarized  In  Section  1.  The 
main  result  is  negative  in  character,  and  shows  that  no  power  of  a  certain 
operator  is  compact. 

The  reader  will  find  the  necessary  background  material  in  Reference  7. 

This  includes  the  definition  and  some  properties  of  Banach  spaces,  the  inte¬ 
gration  of  Banach  space  valued  functions,  and  the  theory  of  compact  operators. 
For  a  more  thorough  treatment  of  the  Fredholm  Theory,  the  reader  is  advised  to 
consult  Reference  8. 
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1. 


The  integral  equation  to  be  examined  is. 


3 '  ",j 


-i 


( i  'A  Y)  '■ 


J(jf 


where 


■■  i.-,  j 

y  "  f y- -y •  -  y>  J 


r  0,  S  is  on  the  Cut-  S  *  and  the  argument  of  the  square  root  is 
chosen  to  be  %  ,  Oi  0<  2  'T ,  0  ■  arg  (s  -  “/4  V  ).  Further,  the  integral 
is  taken  over  all  of  three  space,  and  we  also  assume  that  S,  ^  are  chosen 

»  ,<v  *•» 

so  that  S  *  -  />  t  lc  . 

'  -f 

For  convenience,  we  set 

r  (  \  x  i )  2  r(*\  X)  =  ('xix‘inr  hi  %  -  rj 
>(x,  '  >)  =  /V’;;  -  /i  -  ^Vj^/vVj4 
5  - -  A'*  y+j  -  j-  . 


(1.1)  can  now  be  vritten, 

_t _  <■  —  f  $ t 11  ?  i1-2* 

TTTj  71  ^  ■ 


will  denote  Euclidean  three-space,  while  C, (  tf‘)  will  denote  the 


space  of  continuous  bounded  complex  valued  functions,  having  continuous  bounded 


first  partial  derivatives,  which  are  defined  on  ft*  .  Also  C,  (ft})  will  denote 
the  subspace  of  which  consists  of  functions  which  vanish  at  ,  and 


which  have  vanishing  first-partial  derivatives  at  '■*->  .  For  F  <c  C,  '  ,  ,  we 


define  |/f/|  -  max  j  sup  |  F  J  ,  sup  S  |  ,  sup  fjj, 


i  J  p  ■  1 

8UP  Jv-  ) 


/• 

It  can  be  shown  that  i  1 .  is  a  norm  and  that  with  this  norm,  C,  i J  becomes 
a  Banach  space,  and  also  that  C,  ("0  becomes  a  closed  subspace,  and  hence,  itself 
a  Banach  space.  If  B  is  a  Banach  space,  we  will  use  |_tfj  to  denote  the  Banach 
space  of  bounded  linear  operators  from  B  into  B. 

From  Appendix  I,  it  i3  seen  that,  when  ^  C  T\  j  <  y  "■  • 

(  S'  '4  7 

However,  (p-  )  may  vanish.  Thus ,  when  S  <.  0,  by  ,  ,  ,  7~  , , 

J  T(  y.\) 

we  will  mean  f  *  -  *  where  either  cT  - 

or  else  T*  1  '*  *  .  U»e  limit  in  general  depends  upon  the  side  of  the  axis  from 
which  j-  approaches  S,  but  since  in  both  cases  the  analysis  is  the  same,  ve  will 
not  usually  distinguish  between  the  two.  For  S  <  0,  let  denote  the  integral 


(  y  J  (/  v 
— - - 

«T*(  y  ^ J  ^ 

For  i  *  ,  we  will  thus  understand  equation  (1.2)  to  mean 


'  -  ±  K  >  -  -  , 

h  1  J  St  \,7,) 
if'/)  ' 


where  K  *  lin  •'  . 

»  »/ 
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In  Section  2,  ve  will  show  that  this  is  meaningful.  In  fact,  ve  will 

shew  that  in  both  “d  jCj’^Oj  lim  k'r  exists,  and  is  cospaet. 

(  ■*  o 

This  will  prove  that  (1.2)  can  be  viewed  as  a  functional  equation  in  both  C,  [  ' 

and  C ,  ( i\ J ) . 


In  Section  3,  ve  show  that  when  ,  and  (1.2)  is  written  in  operator 


form  as. 


(  I  -  H  j  ii  -  |  1 

U  V. ) 


that  again,  the  equation  can  be  viewed  as  a  functional  equation  in  both  L,  (  h'*) 
and  t  ,  (•(*)  .  This  time,  however,  ccapactness  is  lost;  ve  will  show  that  in 

both  (C|('f4,)j  and  [C  ( <\5)]  no  power  of  H  is  coopact.  This  result  is  disap- 
pointing;  for  if  some  power  of  H  were  compact,  ve  could  apply  the  generalised 
Fredholm  Theoxy  to  gain  further  information  about  the  solutions  (if  any)  to  (1.2). 
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2. 


r  /  -  ,  \ 

/  n>  \  (  V  f  s')  7 

‘V  >  '  *  *  ^  ) 


In  this  section,  some  theorems  will  be  valid  for  Ct  (  n  '  j  f  ’  C ‘  {  \  1 )  ^ 
as  well  as  for  ^ 


In  most  cases,  the  proofs  will  be 
identical  (except  for  an  interchange  of  syntools)  and  will  only  be  given  for 
*  i ;  *  )  *  however,  this  is  not  the  case,  then  the  additional  reasoning 

will  appear,  ftiroutfiout  the  section,  ^  will  be  negative. 


Define  the  operator  CT(^)  by  w- 7* -  If*  ~  — 

j  ”  *  ✓  i  v-  l.  \  .  /  »| 

i  ;  vf  i  Jr  J  * 


.on 


Thus  for  each  fixed  >,  Cr(^)  takes  the  function  into  the  functii 

(°f  >  }  wV’oT'T”*  *  Clear3y  W'-O  $•  into  <  ,°(n‘) 

Since  vl.c  j}  is  a  linear  operator,  to  show  that  it  is  in  both  *  f- , '  •'  j  J  and 


|l  ,  ( *•  •  we  need  only  shw  it’s  bounded.  An  operator  1-  from  G  into 

is  contained  in  {Oj  if  it’s  linear  and  if  U-'  !  •  sup  j(F{x)  !i  c  ->u  . 


lAmmZ.li  C  ,  ll)  1«  contained  in  both  [C.U1,1]  and  [O  ‘'"M 
PFs  W«  need  only  show  M  O.r  >'i!  '• .  if  y  c  C. ( »v \)  and 
then  ji'  ,  II  .  ,  . _ V.''\  {'ell  1  ,  :l 


!|  *  ' 


Thus,  t.f  },  u  ' 

'j  »  ■*  »  rl>)  •  thu*»  defines  a  function  from  A*’  into  a  Banach  space 
of  linear  operators.  Now,  in  a  Banach  space,  a  Cauchy  sequence  has  a  limit  and 


a 
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hence,  under  appropriate  conditiona,  a  Riemann  sum,  like  ^  CvijJdw.  vill 
have  a  limit  in  the  operator  apace.  This  is  to  say  that  we  can  integrate  certain 
operator- valued  functions  with  the  integral  itself  coming  out  as  an  operator. 

More  generally,  let  (-{;>)  be  a  continuous  Banach  space-valued  function 
defined  on  h  .  We  have 

Theorem  2.2:  If  /*  is  a  measurable  subset  of  R  ,  and  if  ^  1  *■*'  v-4!}-  * 

/* 

then  |  exists,  and  is  an  element  of  the  Banach  space. 


Proof:  Serge  Lang  Analysis  II. 


Let  * )  h  t  l\ 1  ,  and  let  K  t  j  \)  be  a  continuous  conplex-valued 
function  defined  on  lV*  K*  .  How,  if  for  each  ^  ,  given  by 

L  •  '}  ■  ,  la  contained  in  (C,*  <f’)j  (  !<-,  (^*ij  J  ,  and 

as  an  operator-valued  function,  is  a  continue  function  of  y  ,  we  have 
Theorem  2.3:  Let  A  be  a  measurable  set  in  RJ  ,  such  that  j  It  *•  **/ 

A 

and  such  that  for  each  fixed  y)  %  j  )  *■  ty  'J wh*n  ^  h  )  •jj 

we  have  then  i  <l)'»  |  •  j  *  *  j.  V  1»  *  f  ,  where  we  have  equality 

of  both  aides  as  functions  of  )  (  -  '  <  )) 

Proof:  fix  y  ,  choose  ^  '  ‘-.(K5)  ,  and  let  be  a  sphere  closed  in  V  , 


such 


that  j  ''  i  <  ?;/.  >7  <  * 
ir.w 


and 


j  j  ''  j  *  <= 


*  i<V 
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where  g  ^  denotes  the  exterior  of  ^  .  Since  C{y)  is  a  continuous  function 


of  -y  ,  it’s  uniformly  continuous  on  ^  ;  hence,  we  can  partition  ^  into  meas- 

/ 1 

urable  disjoint  subsets  V  (-1  ...nt  such  that  *,  '■i  t  inplies 

//c(x,  -  ch)//  <  £ - —  »  “d  also  80  toat 

/€r 

**  whenever  y  wc  ^  . 

In  each  V,  *  choose  ^  *  ,  and  define  the  step  naps , 


(_C,(R- 


yj  and 


H  :  ^  — t  (J.  (the  coup  lex  nunbers) 


l  'f 


given  by 

Ti 5)  .  c(^;) 

h  v;  -  -‘r.^yry  ) 

We,  thus,  have  for  ^  c  V  that  both  j  (y )  j  * 

and  H  ((>})  ~ 


e 


V'l  I  V 


If  G  *  [•'.  ,(  I\')j  6  |],  will  *not«  (-y  evaluated  at 


For  convenience,  let 


We  have. 


|  G(?)  ^  } 


C-  and  £  <*  j 

V  i  ^ 


C  «,  J^“  (m-,  J-j  0  )  V  ^Uyj  P*J'/ 

/« 

jj  c  -  ni  ii'fii  f  )  p *>  i  7  “  K  <  ^  "J 1  ^ J ^  T  | 
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I,  til/'ll  (  ciDijii  *  il  ,f  c<»-rh)  rfjll)  -j ih'.'iv&il 

J  '  r  C  n  a  -1/1  X 


CfflA  V-^ 

L  =  I  ^  I- 


l  li  ci  '/  6  ■+-  ■—  . 


Mow , 


I-  I 


■t  I.  (  !:l'j  ■*■  |  j  ''h  Vyh)  -1}  ~  j  k ly,  *>j  .j(j)  4r  j  ~ 

r  )/i  /i 


r  ; 

j  (  £  Li  l'  )^<(  ■<  T/ijJ  1  /  -  j  ir^.'OjnWj  I  +-  j/  .) 


VIA 


'■  K  t  ' : 


i 

:,j  'j  i  :i ')  ■  ‘ «  (/■//»)  -  |  K  {  y  li)j  l  ;j  )Ajj  4-  |  (  K  <  -I  *o  ,i  1 ’  '  .)  I 


/•9  A 


Cwl/i 


1  *  '  J  ’v  J  'J-  “1  K  I 


~i  1  a 


1 


j  i  A 


'r  i)-W)  'j  +•  j  n  ( y  '•>)  (jf  ( >Mj  j  ^  ^  t 
i’v.)  A 


Therefore,  |  L-  ;j  L  “  i  K  (  y.  li )  y  <  7 )  ^  |  *  2  11  <1  11  C  *  1  * 

'  ''/I 

But,  since  t  is  arbitrary,  and  is  arbitrary,  we  have. 


Let  |\  i y  li)  and  C(  Jj  be  the  same  as  in  the  previous  theorem. 


We  have  for  a  measurable  set  A, 


Theorem 


os  a 


2.1<  If  |  H  l( 'J)H  *  y  »  then  j  C  (  y)  A  y  is  compact 

/I 

member  of  [ ^C,(«v‘)/  ( £',C( 
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,  3uch  that 


<L  \ 


Proof:  For  fixed  i  ,  is  compact.  For,  let  jl  -f  *  ('h)l!  i'  ,v! 

i* 

then  C (]})  j.  „  -  "f-l  ,  and  there  exists  a  sequent 

/  »,  {'**)  —  •>  .  «*  «*».  li  M  3-Mf  •',!  .fj  -  '■<  *(;/.  •')  II  = 

|  t  -=j  ’  >)  '  II  M?  Mil  —  C. 

Now,  with  the  notation  the  sane  as  in  the  previous  theorem, 

we  have  1 1  (  <•*  <  'J  )  c!  \  ^  T”  ( 5} )  >*  s  j  j  -  (  !  i  ^  ( 3 )  "  :  ■  3  >  *  °  'i 

A  J  nA  ZslA 

\  r«»H»  -i 

si  A  c  ' 

(Which  is  the  stun  of  compact  operators),  it  is  compact.  But,  hence,  [c(  !])-N 

A  7 

is  the  limit  of  compact  operators,  and  hence,  itself  compact. 

Definition  2. it:  Let  be  a  Sanach  space-valued  function  defined  on  the 
interval  [A ,  h  j  ,  and  let  C<=  KJ  .  We  say  that  -f  is  pointwise 
Holder  with  index  at  C  if  there  exists  c  >c  and  /[  >c  ,  such  that 

-L  ( 

for  xc  La>J  IA-t/'-c  implies  jl  {  (*)  ~  (  O  II  ~  ^  I  ^  ^  - 

Theorem  2.5:  Let  ^  ^.v\  be  a  continuous  corplex  Banach  space-valued 


function  defined  on  the  interval  L  *.  '>  j  ,  C  fc  (LMl>)  ,  and  A  either  T  +■  c  <= 

^  ,  £  -  i  <_  ,  then  if  -j  is  Holder  at  V  of  index  G<  M  *  I  , 

l 

I  i'M  {  t ^  *  exists.  Further,  if  f  is  an  end  point  and  if 

(1  «-  ‘  f.  A  -  X 

j{t)  -  C  »  then  again,  the  limit  exists. 
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Proof:  First  we  assume  Cfc  (a*  *»)  ,  Choose  an  interval  X  ^  (°A)  , 
such  that  C  C-  I  (hut  not  an  end  point),  and  li  -f  (A)  *  -ft  It  ^  M  *  ^  ^  *  ~ 


for  K  i  .  Then, 


(  u  +  -  (  1_MJ*  »  (  1  ^  > 
J  -  ~'X  1  A  A  J  >  -  A 


d  *  . 


Nov 


,  r(ii2  -  i:  .  { 


MlA-ci  '.jl'iui* 


1 


I.  c  '  C 


M  is  a  positive  constant.  It  follows  that. 


i  (  ±[‘±iA  -  (  xLxl  j 

,  C-  J  _  A  A  J  _  X  -  l 

^  U  G  x 

Looking  at  j  ~f  1  A  ■  '■  A 


a  . 


a  -  A 


,  we  see  that. 


(  -((*)*<*  __  ,  +(  . 

I  ’  V-  A  •{.  A  -  A  '  -V  x  " 

»ov,  on  r  II  -f±i.:iL‘> I!  *11  ir;.f".!il  »  =  /i 

r  —  X  X  l 


-  A 

(  <  ,  *•’ 

Further,  j  ^  !  *  “  --  I  '• 

I 


.  <-.oo  i  >')-  +  (  0 
A  -  A 


tends  to 


almost  everywhere  on  X  .  Thus,  by  the  dominated  convergence  theorem. 


d  A  r  (  d* 


|  MM  (  I.: OA 

<•  .i.  i  *  -> 


X  -  L 


Let  Is  L>  - 


We  now  need  only  shew  that  |.*m  (  LcJlJ*  exists . 

A  c  J.  a  -  A 

-i 

A  i  J  A  *" 


/  ° 

J  I  _  i  ±>. 

*  - , 


<_  M  l .)*<  t  * 


u  fJ = 

J,  v a  •  rr  '■  -- 
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1  l‘  I  r  l  A  I«  nf'izJL)^  Ar<T„„(  i.;.S 

A  OLlc-O'  U'J 


„  .  A  1  !  ,  -  T  /:+•(_'  c.  S  L 

This  tends  to  log  i  i  ; - -  >  /  - 

l-  U  '*  v. 


This  concludes  the  proof  when  Z  is  not  an  end  point.  In  the 

end  point  case,  chocse  i_  <-  Ll‘-  *• -!  with  T  as  end  point  of  i.  .  The 

analysis  proceeds  as  in  the  first  case,  except  that  since  -f  ('<■)  ~  ^  , 

it  is  clear  that  *  "'l  f  (  ' )  \  •- —  .  exists. 

r.  ••  l 

Definition  2.6:  Let  v.-  «  ’  ,  and  let  -j  ( >}  be  a  Banach  space-valued 

function  defined  on  .  We  say  that  j  (*)  is  uniformly  Holder  of  index 

i>  on  l  if  there  exists  constants  A  ,  i  independent  of  X  in  '■  , 

f.  I;  -f  <*)-  f  1 .1 )  1 1  -  A  1  K'  y'i  ■ 


such  that  > 


<  ,  j  <  l-  ,  and  j  x  -  H  |  <■ 


•  M 


Let  /\  be  an  annular  region  in  (?  ,  given  by  C  ^  &  ~  y 


where  b  >0  .  For  C  real  and  U  V  L  ,  let  equal  either  C  *  1  ^ 


or  r-l.fc  .  If  -fiy)  is  a  complex  Banach  space-valued  function  defined  on 


A  ,  we  have 


Theorem  2.7:  If  L  V  u  ,  if  j.(*)  is  uniformly  Holder  on  /\  ,  and  if  for 


some 


t  >C  1  1 1  II  J '»  ,  then  lim 

A  y~-  AA  L  ->  c-  a 


f  I  1  )  J  U 


exists 


Further,  if  Z  ~  m.  r  r-  ,  then  again  the  limit  exists. 


o  . 
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Proof:  First  assume 


V  <t  L  ^  *>\J  .  For  all  s.  c  A  l1-1-*) 


1  h) 

r-t 


and  further. 


S  ll  till  il  Jh  t  f  II  iilL  | U.  ~ 


\  IliML  j  V  **  1  j*  ^  m  [  HHP! L  J  1 

/  /  AM  |  ^  i  /  ?*-  AN 


I  w.  *■-  A*  I 

where  i'  I  is  abound  c  o  |_£  in  /l  .  So,  by  the  dominated  convergence 

l^-TI 


theorem 


,  i,»i  (  jiiLii  =  f  iilili 

‘  J,°  i  v  x  i  7‘- 1 


Next,  assume 


U  L  LCi  '  ’  ,1  -•/  ;  *  0  >  l?  .  Break  4  up  into  two  regions. 


,4  ,  given  by  cA ^  “  l-  T  t  C~ <■  i>  ;  A  L  given  by  C  u'  f  ^  4,11 


By  the  above,  we  know  that 


i  <  » )  **  3- 


I  mm  (  t!± 

<  *Ak  f' 


C  *  c  x  n 


A 


exists. 


where  we  have  3et 


The  Jacobian  of  this  transformation  is 


1 C  yVl  j  j  -ih)  $""t' 

"a  . »-:T~ 

.j  ^  r’4  (coft  i-My  ^*,c  Ct  H  4  ) 


^  4  J 


/  c  -f  ^  r»  J  ^  y 

Thus,  by  Theorem  2.5,  ///W  •  -- - 5  will  exist  if  we  can  show 

<:  ■-*  c  ja,  y '  -  a 

,  j  n  .  <7 

that  |-  '  (^  I'r)  2  |Vl  (  {  is  a  continuous  function 


I.  / 

‘M. 
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of  r  ,  which  is  pointwise  Holder  at  T  =  "T  .  Let  I"",  be  arbitrary  in  l_&Lj  X  -+<L  J 


Set  ^  a  .J.i  ^  J  m  6  i  i*>  ^  )  Ci-id 

4,=  r,y*  (Ctf6Sm^I<«C  S-  i/cjfjaen, 

|  | '  ± (3)  -/fSfO  f  ...j,  - 

111 'r<  d<.  do 


Now  there  exists  4'  ,  J">G>  such  that  j  cT  and  X  A(£  /I 

4  ,  A-i 

imply  that  j|  -j  (  x)  -  -f  (*>)  II  ~  ^  I  *  ~  *  1  j 

Choose  |  r  '  r ,  j  ^  Cj  .  We  have  then  C  >  I  r '  X  |  ^  -  /  X  —  T,  ^/  = 

1  r**/  independent  of  6  and  (^  .  Thus,  for  |  ■r-  17  J^cA 

j  f  "  11  jh)  -j( y.)//  </  4 J 0  f  7  *4*  /  r  ,/l--  r- 1/1 1  *  _ 

J0  \  "7  Jo  J.,  -I" 

TrVl  ^’A-  r/‘|%  rrv/f  i  r- n  j 


*A 


Ibis  proves  that  <v  A)  is  continuous  on  j_rt  ,  T  M  _i  and  pointwise  Holder 
at  L  .  The  inequality  (r  -  ^/lj  ^ 1 1-  IT,  j  shows  that  f  *  \s  pointwise 
Holder  in  Ca1'/  t  •  It  is  easy  to  show  that  if  ‘jt'Jis  complex¬ 

valued  and  ~f(r)  is  3anach  space-valued,  and  if  both  are  pointwise  Holder 
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at  r,  »  then  ^(r)-j£(r)  is  also  pointvise  Holder.  It  follows  then  that 
r  (vH  is  a  continuous  function  of  lr  which  is  pointwise  Holder  at  r  -  X  . 
Finally,  if  Q.1-  Z  r  o  f  we  proceed  as  in  the  preceding  case;  the 


last  step  being  to  show 


f 


"  r'^Q(r)  cir 
r  -  A 


attains  a  limit  as  £ 


This  limit  exists  by  Theorem  2.5,  because  t  evaluated  at  zero  is  zero 


This  concludes  the  proof. 

Let  K  (y  .ij  be  defined  on  If  *  ^  .  Assume  0(  6  C ,  (  •?  )  ^  C,(i\  )  J  ^ 

where  [_  C(^)j  -/  -  We  have 

Theorem  2.8  If  and  all  its  first  and  second  partial  derivatives 

are  uniformly  bounded  and  continuous  on  n  A  \  ,  then  C(%)  is  a  uniformly 
Holder  function  from  into  ^[C  *  (  (?J)  J  ^ 

j 

Lemma  2.9:  Let  U(ty  be  a  complex  -valued  function  defined  on  If  .  Assume 


the 


partial  derivatives  exist  and  satisfy  |  j  ■*  H  c-',1,!  uniformly 


3  i 

throughout  ^  ,  then  |  *  CM  l  X  ~  I 


Proof: 


I  u(*)|  5  I  -  |  -  w  ( i  •-')  >.  *<>  J  + 

Urhuj  e|UK»-*)-^W.>l+ 
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I  3 1  V.VJO  -  -  I* 


This  completes  the  proof,  since  the  other  two  terms  in  above  inequality  are 
also  less  than  2.  ^  |  A  ~  ^  ! 

Proof  of  Theorem  2.8:  j|  (  OL  }  )  -  ^  ||  -  1 1  ^  ’  <5  ^  ^  ^ 

II  +  'I  t**)  !l  ~ 

,  II  .  |a/u\—  3(<)|  /I  *'  ^0  1 1  ~ 

I  .^h)|  ll  teffr-i)  -  *’(*'*>11  +  1  ^ 

|  p  J)|  II  K(V1)  +  ^  I  3(3)-^X)i  ~ 

„  j  ||  II  K(  pi) -*(*>*’>  H  H  GM  "5"  '*'*1 

where  ^  is  a  bound  on  I'C  and  its  first  and  second  partials.  Now 

I  k  (  -  If  (  *,  *t)  I  4  <>  M  I  *  "  £  I  u”  J 

I  S  {*>^i  ~  ThuJ' 

II  K  (3,*?)  -  l<  L  II  ~  6  *  I  *  "  *  *  * 

Hence  ,|  Gdl-CCOll  «  S'“?  I l(“U)-^)j/l  « 

I  2  01  |  A-  2  I 

This  completes  the  proof. 

We  now  show  that  K  in  (1*3)  is  compact  in  both  [C,  (  (?')  J  and 
.  In  what  follows  K»(flV]can  simply  be  replaced  by 
In  Lemma  2.1,  we  saw  that  •  Taking 
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~  n>efr«  ^  (2.  S')  •  Otpli^i  rl«T  Qfl’g)  'i  uniForm^ 

Holder,  and  hence,  a  continuous  function  of  ^  .  Below,  we  shall  show 
that  j  ||  Ot(^)  II  Jy  <  oo  .  Thus,  from 
Theorem  (2.3),  we  see  that  !<V  =  Jcr($j)cl^ 

Theorem  (2.1*)  then  shows  us  that  jC^-  is  compact. 


Let 


«»)» 


3^j) 


Theorem  (2.8)  shows  that  Oi'jJ  is  uniformly  Holder.  Setting  X  =  Yjfi*%r 
r  -  we  have  KJj. '  (oa-l^)^  *  ( -°f  -3L.  -  j-&- 

In  Theorem  (2.7)  take  *\  *  lv  .  It  then  follows  that  Hr  exists. 

€  -*o 

Thus,  K  is  tM  limit  of  contact  qperators,  and  hence,  is  itself  compact. 

It  remains  only  then  to  shew  that  J  ii  o„(a)ii  *3  <■•  j°. 

n.~  U  ^ II  II  /V' •tr'1)  ~ 


ii  a^'ii 


Ct  ■  Vs)(- 

II  II 


From  Appendix  I  we  see  that  for  some 

I  _ _ _  /  *  — - 

|  •  V^-s|  h( 


A';  a  >  c 

<  _i _ 

M  (  >V4  *  )  ’ 
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Further, 


I  j'i  (  :j'  +  7  *• 


M  1  (  1 


i  '  V/  >',/■!  *  _  ^  K 

777^  rV'-TuV  *4~<V- ,,N<) 

|  /  1.  .V.  j 

where  is  a  bound  on  1  ‘  ■’ '  1 

<  \f  -  4  J  y 

Thus,  there  exists  positive  constants  such  that 


n  \  <  i  ■  i 

i  1  \}  V  4  v- 6 


1  /  -  1  H  <  •>  ) 


IlMftl  1-  ^ - —  Hence,  ||<ri 

( :pdJ  v^-n 


OH  r  /;  ’  , 

f  ii  f  ii‘  • 


-  <-  <.  . 


J  .  C ,  >  c  „ 
•  /  * 


(‘J>  <  J)(  w4*) 

Therefore,  j  1 1  C  r  I  ^  )  IN  y  <•'  *c-'  - 
3,  We  now  write  (1.2)  in  the  form, 

■W  M  -  -  (  H12121-L  «  1LN 

■  Tt  '  )  J  (  hj  t  (  y>')  s  f  V  v, ) 

Let  j“  be  either  >  +  (  c  or  i  *  U  ,  i  *  ^  ,  and  define  Hj  y  • 

1  ^  ( h)  ,'|  1  3  \  '  ^ 

'  r<  vi  *1)  4  •  :j  >  ■“) 

Theorem  3.1s  Mr  is  contained  in  both  {_(-,( i(‘,)  j  and 

Proof:  We  know  that  k',  4  ia  continuous  with  continuous  partials;  hence, 


the  same  is  true  for  Further,  there  exists  positive 

wch  th«  I  (Hrj)c.)|  i 

1  ' '  l  I  L  . l 


.  v.-  .  ||  ,~| 
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s  iT  II  Y  (  ft  Cftf  ♦-  k)  % 

t- \  i ,  — - 

Vi 71 1  c 


-  //j,//  4. 


13. 1) 


where  A  is  some  positive  constant. 


|  k  -  —  (  lllLll.  ...  + 

<)  I  rfy,*))  'f  (  ip7) 

ft  (  -  /— ;  jl— . V(.,  - 

■  lJ!>.  Wj  <  i  )  <  y  >;  f  (y.r)'  ./ 

N 

fr  (p:0Pi)  +-  '  *  % ,  f  -  p  *•  >)  ?P)  ^ 

‘  J  l  rY ;}  1  *■  y,  (  v,  r  ) 

The  differentiation  under  the  integral  sign  could  he  justified,  but 

for  the  sake  of  brevity,  we  omit  the  details. 


Let  i},  be  a  bound  on  |  j  .  The  first  term  in  the  sum 


is  equal  to 


77  1±.±-  -  Ms)  , 

'*  P<’h)  ?('>)  ’ 


and  by 


(3.1),  it  is  bounded  by 


b,  U  jll  L 


r<  M  I 


(3.2) 


The  second  term  is  bounded  by  Hfjj  I  M  *  m)  (  (  -  »  r  f  11 1 )  «* } 

'  J  'V «  V- ~ 

M  jll  (  M  !>|  1  k)  f  - 1  '>1  i  j-  i  L _  i _ 4  .1-  i 

- . - 1  “7  "  y  .  •- . .  \  /  *  *  It  *  '  )  f  y 1  *  l  7 


,  «)a'  (r-‘)'4  ft  • 


1 1  j  <1  (  *■  'ft*  •  '•)  t  (  1  ^ 

f*.  . . ;  (  y  ♦  c ] ,/fc  )(-/'■'  V  •  * )( 't'  * v* ) 

where  is  a  bound  on  *  I  '*•  I  ♦  tyi 


(3.3) 


(rV’0/i 
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(3.3)  is  bound  by  H  ^jll  (  &  I  I  ♦  *0  ^  *-  2.  TTL  (3.4) 

^M’'(YfO/l  i'}1'  ~ 

* 

(3.4)  and  (3.2)  imply  chat  there  exists  a  positive  constant  *v I  ,  such 


that 


j  A.  (  n,j) (h)j  *  AlllAL 


<xr  />■>■! 


’  I  '!,}• 

c 


(3.5) 


(3.1)  and  (3.5)  imply  that  |-j  i-  v-.0\  jj  .  (3.5)  also  shows  that  when 


vanishes  at  oc  .  So,  to  show  that 


H  j  f  s ,  f  r  J  j  it 


remains  to  show  that  when  u  vanishes  at 


1 


h;  m  vanishes  at  >w  . 


For  M 


J 


y*  (.  («) 


write 


r  (  j  v;  7. '  j 


'  T(  I.yt-Cf.r) 

tr 

where  $  is  a  sphere,  centered  at  the 


origin,  in  whose  complement  j  ^  ^ ) j  *  (■  As  before,  we  have  positive 
constants,  such  that 

r  (  S)  j  t  ‘J  )  |  ^  il  '|ij(u  I  *H  ♦  *».) 

1 (  J  ) f  I 


| !  -*J  1 1  (  *  *  »  M  *  *•  > 

M  i* » ,  j  {  *j  * v-  y 


it  |  i)  .. i  •} f  ♦  ! 
r*t «'  i .  J 

Also 


.1  V-')'1. 'Y  ) 


*  < 


<  ‘‘Z  |  <  »)  '  vf  Si 


I  •  $(j  ‘)  *t  5  * 


«  i 


I  (’  ‘  ‘  Hill  I  {  hill  ill 

'  J  N  t»3  <  j  )  I  r,  }  ,)Uy  ’  >  *  ’i 


1 


c.  >  O  ’  i  «  )  . 


•*  K 
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As  a  consequence, 


So  that 


|  (i|r‘i)(''0  |  -  lUVl  +  £'  L- 

given  C  ,  pick  C-  such  that  6  L  *  °/^  ,  and  then  choose  Y  so 

that  |  '//  '  )'  implies  j  -L  c-  ( '')/  *  .  It  follows  that  for  j  ]1  i  >  Y 

This  concludes  the  proof. 

To  show  that  I  •  Hr  exists,  we  first  write  Hr*  ^  ft  r  H  <r  7.  , 

6 

where  H.r:  *  (  1  ^  ,  and  $  is  a  closed  sphere 

centered  at  the  origin,  which  contains  all  zeros  (if  any)  of  ^  (  yi$)  as 

interior  points.  We  now  show  that  H,  ■  I  M  (  exists  and  is  compact  in 

6  -*>  C 

both  [^lOvOj  and  jC  ,  0?0j  .  In  the  following  argument  [C  ,  ( R1  )J  can  be 
replaced  by  |^C  *'  (  >V)  J 

Let  Cr  (5)  y.  •  I  (  '}  i  H  J 

It  is  easy  to  see  that  Crl  t  .  Theorem  2,8  can  be  applied 

0' 

to  show  that  Cj.  (  ^  j  is  a  continuous  function  of  ^  ,  and  hence,  since  ^ 

is  bounded,  j  II  Cr  1  •/)  ||  A  y  ,  Now  theorem  (2,3)  can  be  applied  to 
i  * 

give  ^  C  rW)  ^  Theorera  (2. A)  shows  that  |(rt  is  compact. 

S'* 


Let  CA’g)^-  ■  !  )  .  Tiieorera  (2,8)  shows  that  C  (  i )  is 

r(  ‘j.'jJ 
•>  * 

uniformly  Holder  on  !\  .  Setting  A  *  %  r>  r  ^  ^  and  T  *  Vj  *  ^  * 


we  have  HJ-| 


.  f  cr(jlJ,  ,  f  iMiL ,  1/  (  JiiiL'a-  :%\ 

^  riy.T)  yl-  <  -T)  \  \f~  A 


In  Theorem  (2,7),  Cake  /{  •  .  It  follows  that  |  H./-,-  H(  exists, 


t 


and  is  compact. 


Define 


(  H  ;  a  ))(-'■')  -  (Hz  ^  )(•*»)  “ 


(  r(^  i(L*}  ^ 

5  "  )(  *  U'l.y  tHl'H 

£  o 

,  f  T(  j0(^-s;  ;j  (  *H> 

^  -)  (7..  i,l 


c:\,  iA.V/J  +f.p')  r  f^r) 


'  |  /l^ll  dy _ 

iM  K  ’ .  y  t  c.)  m,  ( 

i  ***  « 


(where  i 1  J  !  :  rre  positive  constants). 

|  J  “  s  I  ‘ijil  (  "  1  ;  <  '  ^ J  *1 _ _  — 

«  ' 

(Note  that  H4  and  may  be  chosen  independent  of  J~  near  5^  .) 

I-T-SJ  II^II  (aril  1  I.J  -ITT  /  j  h  ^ 

- - rl)y'U 

|Vl  A,  1^1  tl)  > 


•m  im  ,  n*  ,i  t 


L  ii  ;j  II 


,  where  L  is  some  positive  constant. 


(3.6) 


Also,  |  .) 


d  (  ;.j ( :i ) } 

J  ,  r( a # V  i(3'0  t(y/0 
c  i 
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(  Lang  Analysis  II)  that  is  contained  in  both  and 

[c,U')|  .  and  since  ||  H*  ~  ||  "  ,  H  (  i!<-'  ~ 


,  A  *  .  _  ,»  * 

L  |  j  -  i  .  L  >  o 


,  we  have  I 1  >vi  M  <r  i  ~  '  V  . 
t  —  c 


Thus,  I  1  ^  llj-  -  H,  *  Mi  *  H  is  contained  in  both  fC  (  I?1)  I  and  [c  ,  ( i\  )  / , 
C  -*c  L  J 

and  further,  H,  is  compact  in  both  operator  spaces. 


In  both  1  C  ,  ( »J4J |j  and  |_C  ,  °(Ry)  J  we  have 
Theorem  3.2.  H  is  not  compact. 


Proof:  Let  “  cl 


-Vo 


It  is  easily  seen  that  there  exists  a  con¬ 


stant  J_  x;  ,  such  that  |7  H  -Jf  L.  .  Now,  for  fixed  ^  ,  we  claim 

CM  (Hl  V.)  (  '?)  -  (/W.  (•■*»)  ,  where  we  take  <1  ?  I  .  For,  in  ^  V  , 

I  — ■»  <.  7 

7"  (  hj  ’i  ( *|j_)  £  '  (  ^  l —  ___  .<£  *  f  *7  J  _ 

f(  ( :i  sJ  cl  ( ^ Vi'-' V  (^7  •  which 

(  fM.M,  are  appropriate  positive  constants.) 


is  an  integrable  function  of  ^  over  (/  $  ,  and 

r  *) 

“*  *■'  Hi-  *>) 

It  follows  from  the  dominated  convergence  theorem  that 

I  "*-i  (H,  (  '*»)  “ 


for  all  W  £  ^  V 


Examining  j-|  L  U  ,  we  find,  for  all  ^ 


12-25 


Ci-  i) 


T  ( °)J 


(  J 


cr  s 


(for  appropriate  h,  f‘  iti  ?  J  >c  ) 

‘■y  )  i  y 


^i(  >*<  V^O^j  (V  4 M 

r(-h)  '  d 


f  __jjj _ 

J  .  *?lu  ij  (  / i  J  ‘ ) 


xv  lK  if. 


l-i  i  lvi  3 


r  i  <>i4  cl  (j  ,-1 1-  t3  r 


•A  J  ^  *,*-*<  *)(  » '■' 


where  ^  is  the  radius 


of  >  .  This  last  term  is  equal  to  'f’('h)  1  ^ 


XJ 


«*M.J  p 


r  d  r _ 

(rS  V<< 


r*  r  r 


^  II  1  ^  ’‘0 

i  i ,  i-i  j  U  ^  **  ’-1  v‘^0 ri  (■*■") 

(where  on  L  ^  ^  c  L  ) 


>e-> 


i  fTT(^) 

'3  v  '  '  'p 


XJ 


,  ,  _  -F  IT  Tf  -J 


'Ml'M,  l.  *  <  'J  )«('  ^  V'"'*)  f‘.  T,,(  »  1  .) 


Aii  1  n  »i 


tr  t  (  ■•>.) 


ij/  _  /\  r  i  f  m  ■ )  _1 


V  V  1  ^ 


l  Vi  C 


\  1  1  <■  L. 

K  yo. 


r  1  y  ' 


Now,  *  1  ,V| 


_  M  "i 


( "?i  —  «  IVh,  1  "‘I-*  " 

Therefore,  for  large  j  **!  j 

Kh.j ’)<•>.) I  >•  (fi'ftj 


*■  -i  a 


We  will  have  shown  that  |-|^  ,  and  hence  H  ,  it  not  compact  if  we 
show  no  subsequence  j  H*  ^«<j  converges.  So,  assume  Ha  converges' 
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£< 

to  say  ^  .  i.e.  ({  Hx  -  -f  ij  **  O  .  Now,  since  Hz  6  C  t  (  R*)^  jl 

Bust  also  be  in  C,  ,  But,  for  each  ^  we  must  have  ■ 

I  "'M  (Hi  *  (^1-5  )  (^)  .  This  implies  that  for  large  'a7 

I  f(^j  ■  j(^*  ^  )(  Oj  ^  )  ^T"  ,  which  contradicts  the 

fact  that  C  Ct  ( l?JJ  .  Thus,  does  not  converge.  This  concludes 

for  proof. 

In  both  [__C  ,  ( /?sjj  and  [C,  °( R*)J  ,  we  have 
n 

Theorem  3.3.  If  ir?  I  H  is  not  compact. 

Proof:  We  use  the  fact  that  if  A  is  a  compact  operator,  then  for  any 


bounded  operator  If  ,  A  if  and  if  A  are  also  compact.  It  follows  that 

o  n 

|-|  -(Hi  Ml)  will  not  be  compact  if  and  only  if  is  not  compact. 


Let  j  b«  a  positive  integer. 

n  -v« 

Theorem  3.2  we  have  this  for  -  I  .  So 
it  holds  for  j  * 

Ut  fnj  (■’?)  -  ft  ?<$»)(*)) 

(Hx  ^)h)  =  (^*-  ')/\)  " 


We  first  show  that  for  fixed  ^ 
where  ^  )s  I  .  In  the  proof  of 
assuming  it  holds  for  -j - i  ,  we  show 


and  ■  (f-ij.  . 
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kov  iifnr'ii *  ii  nr^H  - « h;"u i.j-n  *  l j 


where  Lj  ,  is  a  positive  constant.  Thus,  for  ^  c-  Cf  5  and  some  M,  Mt 


L  ">  C  we  have  - 


_  -T(h)L« 


/-J) 


which  is  integrable 


over 


€i 


Further,  for  ..ch  y  c  C  i  r(3,<,)+h.O  r(S,'.jf  tjM 

Thus,  by  the  dominated  convergence  theorem,  for  fixed  ^ 

(tiSn i)(-h)  --Iih  (  Hi  -f  H  1-1  )  r  (Hxfjtjh) 


i  'vi 


Secondly,  we  show  that  ('ll  ^  *<t  cvo  ...  .  We  have,  in  Theorem  3.2, 

proved  this  when  j  “  I  . 

Now,  since  >C  in  C  S  »  we  must  have  (/-|z  ^ 

and  similarly,  it  follows  that  (Hi  ^ c  "  »•  know 

from  the  proof  of  theorem  J.2,  that  is  bounded  away  from  G  at  *->  . 

I-1  V 

So  assume  that  Ht  jj  is  bounded  away  from  C  at  x>  ,  It  must  then  be 
bounded  away  from  soro  throughout  R  ;  that  is,  there  must  exist  <*j~i 
such  that  ^  (Ht  ^  )(^J  ~  0 

(Hj^i -ijf-lj  r  *1 V,  >0  .  Ih.r.(.r., 

1  4 

Now,  if  |J  M C  *  th#n  /  would  have  to  be  Ht  ^  »  •nd 


>c- 
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sC  the  same  time  j  t  (,  .  Since  this  is  impossible,  H*.  is  not 

compact. 

This  concludes  the  proof. 
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Appendix  I 


Leona:  Let  A,  Q t  R  and  $  be  real  or  conplex;  then  if  either  l\e  S'  <  C. 

or  I'h  *  tc’ ,  there  exists  <--*-/  M.(  <HA >< ,  such  that 

m,  (  n\c\)  -  |  a1*  /i*  r?  i  n1-  ■»  |  *  M  i  (  A  ^  RV*  a  } 

i. 

Proof:  /\  h  /'  *  <i  *  »i  '>-/(‘(  i]L-  2.  I  4 1  |  IJ I  -  (  I  *  I  -  I  Ul  j  ^ 

C  with  equality  possible 

only  if  j  A  |  ~  1  G  I  .  Now,  if  |A|^|tlj  then  /\l+  /i  ■  i  0  -  jjffnu- 

,  where  £  is  angle  between  A  and  3  , 

On  l\  *  let  A  *  A*  Q+  0  ,  Let  O’  be  the  surface 

of  the  unit  sphere  in  Euclidean  R  *  .  Then  if  ( 4 ,  and 

I  A  *  ■  I  &  I  .we  have  /A  I  m  />_  .end  hence,  /  A  |  *(  *  *  L  *-  S  t  )  f  C  , 

It  follows  that  for  i? J  4  U  F  (A,i))  >C  , 

Tlius,  since  \J  is  conpact,  there  exists  such  that  for 

U  f"  ,  Further,  for  arbitrary  we  have 

(VO  •  (  A  At/  A  where  >1.  t  0.  “  I  XsA,0 
\(  A>*\)  *  A  *  ( A< .  ft.  J  .  Hence,  Am  t  f  (a,  ft)  *  A  M 
£>- 

For  /7»,  0}f  t  ,  we  have  i*l  •*  i~(A,nj  *  M  ,  Now,  choose  a  >  c  , 


we  have 


For  |  (A,R)  j  L  , 


I  ?(*><*) I  ^  f(A,n)  i  s  1 

A1 1  Q“  A'*  il' 


lvl  +  I  S  *  ,  and  if  L  is  large  enough  we  have  F  (  Aj  _  —  I  S 


f-.rto.  ALtQxiu 


/* 


jr(Atu)'i  I  md  |  f ( a ,ii)  -  s’  j  ^  r( a, i?j  - 


/al  *  qx*  a. 


!il  -  11L  . 

i  2.  u 


A1  ^ 


i  (  /\‘  ♦  2  (  /\'i  >y) 


Now  choose  L.  so  that  I'M  +JLli  *2  rH  end  -23-  ~  i-i-i  ^ 

L  2.  1  L.  ' 


Us  thus  have  for  j  (A ,  •"? )  | 


^  I  F(A,tf)-  Si  *  2  »v| 


J  4 


aA 1  t  «3  *  -i  a. 


For 


mi  I  F(A,4)~  1  I  _i  ^ 


/A  X  ♦  Q  '  i  u 


for  som  ivj ;  i n  >  t .  Cat  Mx  *  Pio*  [JlAi,  nj 

rM,  *  r».n|7»»^#  '7i  I 


Va  have,  thus,  for  all  (A  ,  ft ) 

M*‘i  a1* «)  4  •  fM'0)  * 


I  t  M,  (  /*'■>  «J 
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